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INTRODUCTION 


I.  LASER  TECHNOLOGY  AND  PROPAGATION 

For  CW  laser  propagation  through  a  continental  aerosol  distribution,  wp  can  approximate 
the  rms  deflection  ■'fa  laser  beam  relative  to  the  diffraction  spreading  '.id  find  the  effect 
to  be  small  for  a  1-kW  beam  of  2-cm  diameter  =  x  10  Z)  over  a  kilo¬ 

meter  path.  However,  since  it  scales  as  power/diameter,  it  can  be  significant  for  larger 
devices. 

Further  results  have  been  obtained  in  the  laboratory  experiment  on  CW  laser  thermal 
blooming.  The  addition  of  a  nonoptimum  phase  profile  to  the  near-fieid  laser  beam  re¬ 
sulted  in  a  50-percent  improvement  of  the  far-fie’.d  peak  irradiance.  The  peak  irradi- 
ance  on  target  was  observed  to  decrease  when  the  laser  powe”  was  increased  beyond  a 
critical  value.  An  improved  method  for  adding  the  desired  phase  profile  to  the  laser 
beam  involves  a  pair  of  cylindrical  lenses  and  a  pair  of  spherical  ienses. 

The  effect  of  pulse  trains  on  laser  propagation  has  been  studied  for  a  few  representative 
cases,  showing  that  for  the  same  total  energy  thermal  blooming  is  relatively  insensitive 
to  the  pulse  forms. 

The  effect  of  a  precursor  preceding  the  s  .cot.l  pulse  by  more  than  three  hydrodynamic 
times  is  being  studied  using  a  double -pulse  program.  Results  indicate  that  lie  second 
pulse  can  experienct  subs'anMal  blooming  when  the  energy  of  the  precursor  is  large. 

We  have  developed  a  propagation  code  tir.t  includes  one  step  of  time  integration  as  well 
as  three  space  dimensions.  Calculate v.S  of  blooming  for  square  apertures  have  been 
carried  out  and  compared  with  circular  apertures. 

Numerical  computation  has  shown  evidence  of  Bruerkner-type  dynamic  h.st&biiities  in 
thermal  blooming  of  pulsed  laser  beams.  Furthe>  exploration  is  planned  to  determine 
whether  the  effects  of  turbulence  can  oe  amplified  to  produce  these  instabilities. 

Multiple-scattering  calculations  have  been  Carried  out  for  propagation  of  laser  pulse 
through  a  path  in  a  fog  cleaied  by  a  precursor.  The  rms  oeflection  of  this  second  pulse 
is  proportional  to  the  third  power  of  rm,  the  peak  valu  .•  of  the  fog  drop  size.  Thus,  for 
fogs  of  considerable  water  density  nnd  rm  ~  2  pm,  spre  ading  can  be  significant.  The  ex¬ 
isting  Bradley-H jrrmann  propagation  code  is  beir.g  adapt*  ♦-  it  elude  propagation  through 
a  fog  or  aerosol. 

The  absorption  of  a  laser  beam  by  a  saturable  medium  has  bt  hi  examined.  The  steady- 
state  absorption  of  a  focused  beam  by  a  simple  three-level  system  is  shown  to  be  e^uiv  • 
alent  to  that  of  a  saturable  absorber.  This  differential  equation  has  been  solved  miner' - 
cally  for  several  different  conditions.  The  effects  of  focusing  are  significant. 

Saturation  'txperiments  have  been  performed  with  a  focused  10.6-pm  laser  beam  in  a  5-tn 
cell  containing  500  torrs  CO-  and  250  torrs  argon.  For  these  conditions,  the  linear  ab- 

^  .4 

sorption  coefficient  was  1.8  x  10  cm  and  the  transmission  of  laser  energy  as  a  func¬ 
tion  of  incident  intensity  indicated  ‘hi*  saturation  intensity  was  from  1  to  2  x  105  W/cm2. 
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Breakdown  experiments  were  done  at  1.06  jim  arid  indicate  ~  A  '  scaling  for  boih  clean 
and  dirty  air  thresholds. 

Growth  dynamics  for  particle-induced  air  breakdown  were  studied  to  obtain  plasma  growth 
rates  and  absorption  depths  as  a  function  of  intensity.  At  lowrer  intensities  the  velocities 
increase  linearly  proportional  to  the  intensity  1;  at  higher  intensities  the  axial  velocity 
varies  roughly  as  1*^3  in  agreement  with  planar  LSD  wave  theory. 

Measurements  of  peak  surface  pressure  and  impulse  delivered  to  aluminum  and  carbon 
targets  are  reported  in  which  two  CC>2  laser  pulses  separated  by  a  variable  time  delay 
were  used.  Delay  times  for  recovery  of  peak  pressure  and  impulse  to  single-pulse 
values  vere  observed  to  be  about  1  msec  for  the  conditions  of  the  experiment. 

Results  are  presented  from  an  experimental  examination  of  the  plasma  environment  pro¬ 
duced  by  irradiation  of  an  aluminum  target  with  two  10.6-pm  laser  pulses  separated  in 
time.  Plasma  diagnostics  include  interferometric  measurements  of  neutral  and  electron 
density,  cross-beam  10.6--pm  absorption  measurements,  and  fast  streak  photography  of 
the  luminous  front  motion. 

An  empirical  relationship  has  been  determined  of  the  time  development  of  the  radial 
growth  of  the  vapor  sphere  surrounding  a  water  droplet  that  has  been  irradiated  with  a 
nigh-intensity  10.6-pm  laser  pulse. 

For  a  cylindrically  symmetric  density  variation  in  a  material  of  uniform  composition,  a 
numerical  approach  based  on  linear  splines  is  used  to  obtain  the  density  profile  from  in¬ 
terferometric  data.  The  accuracy  of  this  method  has  been  examined  for  various  test  cases. 

The  dual  500- J  e-beam  laser  is  in  the  final  stages  of  checkout.  The  device  is  complete 
except  for  the  automation  of  certain  auxiliary  sequencers.  Each  laser  has  been  tested 
independently  and  has  operated  successfully  with  good  beam  quality.  Double-pulse  oper¬ 
ation  hi  j  been  achieved.  The  two-pulse  laser  facility  has  been  constructed  to  rtudy  prop¬ 
agation  and  target  interaction  phenomena  for  an  arbitrary  interpulse  separation. 

Calculations  indicating  the  feasibility  of  an  electric  discharge  H2~HF  laser  are  described. 

Components  have  been  ordered  for  a  60-J  e-beam-initiated  HF  laser.  An  intense  elec¬ 
tron  beam  \?'ll  be  propagated  nearly  col'Inearly  with  the  laser  axis  through  a  mixture  of 
and  C2H^.  The  e-beam  will  be  confined  byan  axial  magnetic  field.  Computer  stndies 
are  being  made  to  aid  in  obtaining  maximum  uniformity  of  this  field.  The  laser  is  ex¬ 
pected  to  be  operational  next  reporting  period. 

Performance  of  our  transverse-discharge  HF  laser  was  improved  from  0.13  J  to  2.25  J 
output  by  (1)  minimizing  the  combined  inductance  of  the  Marx  band  and  laser,  and  (2) 
changing  the  hydrogen  donor  from  Hj,  to  C2H^.  Further  increases  can  be  expected  by 
increasing  the  input  electrical  energy  with  a  commercially  obtained  Marx  bank  tc  be  de¬ 
livered  shortly. 

Parts  have  been  ordered  for  construction  of  a  10-J  H2-F-,  flashlamp- initiated  laser. 
Assembly  will  begin  shortly. 

An  HF  kinetics  program  is  being  developed  for  the  purpose  of  predicting  the  performance 
of  our  high-energy  HF  lasers.  The  Runge-Kutta  method  is  being  used. 


via 
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Detailed  results  on  "C  02,  C  O 2<  and  C  02  isotope  lasers  will  be  published  in 

the  February  1974  issue  of  the  Journal  of  Molecular  Spectroscopy.  Experiments  with 
12  16  18  13  16  18 

C  O  O  and  C  O  O  isotopes  and  improved  grating-cor.trolled  lasers  will  be 
continued. 

II.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 

A  descripUon  of  the  operation  of  the  new  engineered  model  of  the  fast  scanning  MTF 
interferometer  is  described.  A  technique  is  outlined  to  use  this  interferometer  both  for 
accurate  angle-of-arrival  measurements  and  as  a  new  type  of  nonim'.ging  tracker. 
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I.  LASER  TECHNOLOGY  AND  PROPAGATION 


A.  CW  PROPAGATION 


1 .  CW  Propagrtion  Through  an  Aerosol 

As  a  first  considei  ation  of  the  propagation  of  a  CW  laser  beam  through  an  aerosol  we  take 
the  case  of  multiple  scattering  of  the  radiation  through  a  density  of  index  spheres  whose  size  we 
assume  to  be  sufficiently  large  so  that  the  geometric  optics  approximation  is  appropriate.  The 
mean  deflection  per  collision  is  €  =  An/n  and  the  rms  angle  of  deflection  through  a  path  length  L 
of  index  spheres  of  density  Nr  is  given  by* 

Gr  =  L</Z  [Var2(t)|l/Z£ 


where  aft)  is  the  radius  of  the  vapor  sphere  surrounding  an  aerosol  particle  of  initial  radius  r. 
The  average  squared  deflection  of  the  beam  due  to  a  number  density  Nr  of  particles  of  radius 
between  r  and  r  +  rdr  is  ther. 

yZ(r)  =  4/9  L3Nr*a2(t)  (An/n)2  . 

If  we  assume  that  the  aerosol  particles  have  reachea  steady  state  and  all  of  the  heat  absorbed 
from  the  incident  *"adiation  diffuses  out  into  th'  index  sphere,  ther  the  mean  change  in  index  of 
refraction  of  the  spheres  is  given  by 

ATair  I(TraZer>  t/(4/3  ,arpcn> 

An  =  (n  -  1)  =  in  -  1)  - — — - L__E 

o  o 

where  AT  .  is  the  temperature  change  of  the  air  composing  the  index  sphere,  n  is  the  index 

alr  -4  2 

refraction  of  air  (n  -  1  “  3  x  10  ),  I  is  the  incident  radiation,  *raer  is  the  absorption  cross 

section  of  the  aerosol,  p  the  air  density,  and  Cp  the  heat  capacity  at  constant  pressure  for  r.ir. 

Since  the  diffusion  solution  ie  appropriate,  the  index  sphere  expansion  is  given  by 


a‘(t)  =  Kt 

where  «  is  the  diffusion  coefficient  for  air  (k  =  0.2  cmZ/sec).  Hence,  we  'ind  that  the  rms  de¬ 
flection  is  independent  of  time  if  the  beam  cross  section  is  uniform  in  intensity,  even  witl  a 
wind  present, 

_  vL»<-.> 

n  4«par(t)  y/(y  -  1) 

where  p  is  the  ambient  air  pressure  and  y  is  the  specific  heat  ratio  for  air.  For  the  aerosol 
distribution,  then, 

lV»rr . /.I2  » 

y  (r)  - 2~2 — 7 - 

p  iy/y-1) 


1 


and 


total 


-[C" 

*■  mi 


max  2,  .  , 
y  ( r)  dr 


i/2 
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For  an  oxamole  o*  a  typical  continental  aeroso'  we  can  write  N(r'  =  3r  Nr  where  N  is 

o  min  o 

the  total  particle  density.  This  gives  a  value  of  rms  deflection 


ytotal 


3i/2„  3/2N  /2  „3/2 

- £> - - - 

2*p(y/y  - 1)  n 


,  3  i1/2 

(r  r  .  )  ' 
max  min 


where  r  »  r  .  . 
max  min 


The  diffractive  spread  of  the  same  bean  is  y^f  =  l.x/D,  where  D  is  the  beam  diameter, 
so  the  relative  spread  due  to  the  density  of  the  fog  compared  with  diffraction  is 


ytotal 

ydiff 


3i/2pL1/2Ni/2  r3.  )‘/2 

o  max  min' 


n  -  I 


D*Xp<  -  1)  x1'2 


with  P  being  the  laser  power. 

For  J,  =  i  km,  D  =  2  cm,  N  =  103  cm  3,  P  =  1  kW,  and  X  =  10.6  pm,  r  =  20  pm, 
-  o  max 

r  ar  4  x  10  pm,  p  =  i  atm, 
mm 

ytPtal  =  1.4  x  10-2 
ydiff 


At  this  level  of  power  and  this  beam  size,  the  spread  due  to  aerosol  scattering  is  not  important; 
but  as  CW  lasers  scale  up,  the  effect  could  become  a  significant  source  of  beam  distortion.  For 
nonuniform  inte»._,.ty  cross  sections  and  in  the  presence  of  a  cross  wind,  analysis  indicates  that 
more  severe,  asymmetric  distortions  would  be  effected. 

L.  C.  Bradley 
H.  Granek 


2.  Rec ilts  of  Laboratory  Experiments  on  CW  Thermal  Blooming 

Some  further  results  have  been  obtained  from  the  laboratory  experiment  on  CW  thermal 
blooming.  The  first  set  of  results  was  reported  in  the  previous  Optics  Research  Report.3  Also 

o 

given  in  Ref.  3  is  a  con  olete  description  of  the  experimental  setup.  In  brie!,  a  focused  5145-A 
laser  beam  is  slewed  by  i  rotating  mirror  through  an  absorbing  gas  (NOz  buffered  by  air),  and 
the  amount  of  thermal  blooming  Is  measured  by  having  the  slewed  beam  pass  over  a  pinhole 
array,  behind  which  is  a  detector.  Certain  problems  that  nave  been  encountered  -  thermal  con¬ 
duction.  nonuniform  gas  mixing,  and  laser  multimodirg  -  affect  the  experimental  results  re¬ 
ported  here.  These  problems  are  being  corrected. 

(a)  One  of  the  main  goals  of  the  laboratory  experiment  is  to  verify  that  phase  corrections 
added  to  the  near  field  of  the  laser  beam  can  decrease  the  amoint  of  thermal  blooming  caused 
by  inde?-of-refraction  variations  in  the  absorbing  medium.  A  set  of  zero-power  cylindrical 
lenses  will  be  used  to  add  part  of  the  desired  phase  contours.  At  the  present  time,  a  set  of 
cylindrical  lenses  with  a  good  surface  figure  (X/10)  is  not  available.  Some  daia  were  obtained, 
however,  with  a  ae*  of  lenses  with  -IX  surface  figure. 
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Fig.I-1.  Laboratory  results  which  show  a  decrease  in  thermal  blooming 
when  an  appropriate  phase  profile  is  added  to  the  near-field  laser  beam. 

Np  is  the  distortion  number,  which  is  inversely  proportional  to  the  cross- 
wind.  Ip  is  the  peak  irradiance  in  the  focal  plane. 

The  results  are  shown  in  Fig.  1-1.  Tl.s  peak  irradiance  in  the  focal  plane  was  plotted  as  a 
function  of  distortion  number  for  two  cases:  when  no  phase  correction  war  added  to  the  beam 

and  when  the  phase  correction  in  the  upper  right-hand  comer  of  the  figure  was  added  to  the  beam. 
In  this  experiment,  the  distortion  number  was  changed  by  varying  the  crosswind,  so  that  the  dotted 
line  would  result  in  the  absence  of  thermal  blooming.  At  low  values  of  Np,  where  little  or  no 
thermal  blooming  occurs,  the  phase -corrected  beam  has  a  lower  peak  irradiance  than  the  uncor¬ 
rected  beam  because  the  phase  correcti  -t.t  represent  aberrations  which  prevent  the  beam  from 
focusing  into  a  small  spot.  As  the  distortion  number  is  increased,  however,  the  aberrations  in 
the  absorbing  gas  caused  by  thermal  effects  begin  to  cancel  the  original  aberrations  in  the  phase- 
ro  re :ied  beam,  so  that  for  Np  >  170,  the  phase-corrected  beam  his  r  higher  peak  irradiance 
thar.  the  uncorrected  one 

On  comparison  of  the  experimental  points  to  the  dotted  line,  it  is  evident  that  the  amount  of 
improvement  is  small,  about  1.5  dB  maximum,  whereas  Bra  iey  an<j  Herrmann*  predict  that  as 
much  as  6  dB  improvement  should  be  possible.  The  main  reason  for  the  relatively  sr.  all  amount 
of  improvement  is  that  the  actual  phase  correction  is  far  from  the  optimum  one  calculated  by  the 
propagation  code*  and  shown  in  the  lower  right-hand  corner  of  Fig.  1-1.  More  improvement  is 
expected  once  lenses  with  a  good  surface  figure  are  obtained  and  used  in  future  experiments. 

(b)  The  theory  developed  by  Bradley  and  Herrmann*  predicts  that  the  peak  irradiance  on 
target  will  actually  decrease  as  the  laser  power  is  increased  beyond  ..  certain  critical  power  Pc. 
This  prediction  was  tested  experimentally  in  the  lnooratory,  and  the  results  are  shown  in  Fig.  1-2, 


t  The  thermal  blooming  distortion  number  is  Nd  *  (l/fpcp)  (8t/8T)  (k<*PR/av0),  where  f,  p, 

Cp,  a,  and  T  are  the  dielectric  constant,  density,  specific  heat,  absorption  coefficient,  and 
temperature  of  the  gas,  respectively;  k,  P,  and  a  are  the  wave  number,  power,  and  1/e  power 
radius  of  the  laser  beam,  respectively;  R  ip  the  range;  and  is  the  crosswind  velocity. 


) 
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Fig.  1-2.  Laboratory  results  which  demonstrate  that  the  peak  irradiance 
on  target  can  actually  decrease  when  the  laser  power  is  increased.  The 
discrete  points  were  measured.  The  solid  curves  were  calculated  from 

theory:  I  =  (P/Aq)  e  ^**c. 


y  |u->-mu| 


Fig.  1-3.  Phase  profile  calculated 
to  give  maximum  decrease  in  ther¬ 
mal  blooming  for  a  typical  set  of 
laboratory  parameters.  The  con¬ 
tour  linen  i-  thir  and  succeeding 
figures  are  marked  in  units  of 
wavelength. 


»  la-i-Miil 


Fig.  1-4.  Phase  profile  provided  by  a 
plano-convex  and  plano-concave  cy¬ 
lindrical  lens  pair.  The  radius  of 
curvature  is  14  cm,  the  separation  is 
0,001  cm,  and  the  axes  of  revolution 
are  at  an  ingle  of  0.08*. 
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along  with  the  theoretical  curves.  It  is  to  be  noticed  that  (or  all  four  values  of  vq,  1  does  reach 

a  maximum  va'ue  and  then  decreases  as  P  is  increased.  As  expected,  the  amount  of  blooming 

increases  as  the  crosswind  v  is  decreased.  However,  ai  the  lowest  crosswind  values,  the 

o 

amount  of  blooming  is  independent  of  v  .  This  independence  may  be  attributed  to  the  fact  that 

in  the  experiments  reported  here  conduction  is  the  dominant  mechanism  of  heat  transfer  at  low 

values  of  v  . 

o 

Bradley  and  Herrmann  predict  curves  which  may  be  fitted  by  the  scaling  law 
,  P  -p/pc 

V  =  7T  e 

r  o 

1  is  the  peak  irradiance  in  the  focal  plane,  P  is  the  laser  power,  and  Aq  is  the  beam  area  in 
the  focal  plane  in  the  absence  of  blooming.  is  the  critical  power  at  which  1  reaches  a  maxi¬ 
mum  and  is  given  by 


2v 


P  =  12.4 
c 


:pcp 


akr 


d«73T  ln(l  +  N  ) 


This  theory  does  not  include  the  effects  of  conduction,  which  are  obviously  present  in  the 

experimental  results.  A  nmple  method  of  including  conduction  which  gives  correct  order-of- 

magnitude  results  follows.  In  the  expression  for  P  ,  heat  transfer  by  forced  convection  is  taken 

into  account  by  *he  factor  v  .  To  take  into  account  additional  heat  transfer  bv  conduction,  one 
J  o  - 

can  add  to  v  an  effective  wind  due  to  conduction,  v  =  a/t  where  t ,  is  the  time  required  for 
o  a  2  cl  d 

conduction  to  balance  the  heat  input,  is  equal  to  pCpS  /K,  where  K  is  the  thermal  conductiv¬ 
ity  of  the  gas.  The  curves  in  Fig.  1-2  were  calculated  using  (v  +  vd>  =  v  in  the  e:  press  ion  for 
Pc,  with  vd  =  1  cm/sec. 

Quantitatively,  the  experiment  and  theory  agree  to  well  within  a  factor  cf  two  in  the  regions 
where  the  theory  is  valid  (Pc  £  1.5P-  Ref.  4).  Qualitatively,  however,  the  agreement  is  not 
good.  That  is  to  df.y,  the  s’  ..pe  of  the  experimental  curves  does  not  follow  a  Pe  p/pc  depend¬ 
ence.  The  experimental  curves  peak  sooner  and  fall  off  more  slowly.  The  main  cause  of  the 
discrepancy  is  probably  the  facile  manner  in  which  conduction  was  inc-'rted  into  the  theory, 
especially  in  view  of  the  fret  that  the  discrepancy  appears  to  increase  ?.s  tile  conduction  effect 
increases.  Further  experiments  In  which  conduction  effects  will  be  reduced  are  being  planned. 

(c)  Improved  methods  for  adding  phase  profiles  to  the  laser  beam  have  been  devised.  Fig¬ 
ure  1-3  shows  a  phase  pi~t  computed  by  the  propagation  code  of  Bradley  and  Herrmann.  This 
phase  plot  Is  predicted  to  cr.use  an  optimum  reduction  in  thermal  blooming  for  a  typical  oet  of 
parameters  in  the  laboratory  experiment.  The.  main  feature  of  the  plot  is  a  saddle  point.  A 
saddle  point  can  be  obtained  by  a  plano-convex  and  plano-concave  cylindrical  lens  pair  of  matched 
curvatures  by  rotating  the  axes  of  revolution  slightly  with  respect  to  one  another.  Figure  1-4 
shows  the  phase  contours  that  result  when  the  axes  of  the  lensei  make  an  angle  of  0.08*. 

Figure  1-5  shows  the  difference  between  the  contours  of  Figs.  1-3  and  1-4.  This  difference 
can  be  approximated  by  a  matched  pair  of  spherical  lenses,  one  A  which  is  tilted  out  of  the  plane 
perpendicular  to  the  optical  axis.  Figure  1-6  shows  an  example  of  the  contours  avai'°ble  from 
the  spherical  lenses.  Figure  1-7  shows  the  difference  between  the  contours  of  Figs. , -a  and  1-6. 
The  surface  of  Fig.  1-7  differs  from  flat  by  less  than  A/*0  over  the  whole  central  portion  of  the 
beam,  within  which  more  thar.  70  percent  of  the  energy  lies  (assuming  a  Gaussian  intensity  dis- 
tributlon  truncated  at  the  Ve  power  diameter  by  the  lenseo). 
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Fig.  1-5.  Residual  phase  after  the  phase 
contour  of  Fig.  1-3  has  been  subtracted 
from  that  of  Fig.  1-4. 


Fig.  1-6.  Phase  profile  provided  by  a 
plano-convex  and  plano-concave  spher¬ 
ical  lens  pair.  Radius  of  curvature  is 
25  cm,  the  separation  is  0.2  cm,  and 
one  lens  is  tilted  3*  out  of  the  plane 
perpendicular  to  the  optical  axis. 


Fig.  1-7.  Residual  phase  after  the  pti  ve 
contour  cf  Fig,  1-5  has  been  subtracted 
from  that  of  Fig.  1-6. 


Thus,  a  oair  of  cylindrical  lenses  in  conjunction  with  a  pair  of  spherical  lenses  can  closely 
approximate  the  desired  phase  contour  of  Fig,  1-3.  In  addition,  the  contours  of  Fig.  1*4  show 
almost  pure  astigmatism,  while  those  of  Fig.  I-fc  show  primarily  coma,  so  that  these  two  aberra¬ 
tions  can  be  separa.ely  added  to  the  beam 

D.G  Fouche 

B.  PULSE  PROPAGATION 
i .  Theory 

a.  Effect  of  Pulse  Trains 

(1)  Gas  Dynamics  of  Pulse  Trains 

We  have  performed  calculations  with  the  linearized  hydrodynamic  equations  in  order  to  study 
the  gas  density  change  caused  by  pulse  trains.  These  are  compared  to  the  gas  density  change 
due  to  a  single  pulse  with  the  same  total  energy.  Figure  1-8  shows  the  gas  density  as  a  function 


Fig.  1-8.  Gas  density  as  a  function 
of  radius  for  two  times  for  the  pulse 
train  (o)  and  for  a  single  pulse  (x).  < 


of  radius  for  different  times  for  the  single  pulse  and  the  pulse  train  for  the  following  conditions: 
t2/tj  =  5  ,  Tp/t2  -  5  ,  Tp  =  500  psec  , 

whei\.  tj  is  the  duration  of  a  single  pulse,  t2  the  period,  and  Tp  is  the  duration  of  the  pulse  train. 
The  beam  shape  is  sn  infinite  Gaussian  with  i/e  power  radius  of  6  cm, giving  a  hydrodynamic 
time  tj^  *  200  psec.  The  peak  pulse  energy  density  absorbed  by  the  air  for  each  of  the  two  pulse 
forms  is  8.5  x  lo'2  J/cm3.  As  can  be  seen  from  Fig.  1-8,  the  gas  density  changes  are  not  very 
different  at  the  end  of  the  pulse  train  (t  =  T  )  and  we  expect,  therefore,  that  the  pi  jpagation  is 
insensitive  to  pulse  forms. 
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NORMALIZED  IRRADIANCE 


Fig.  1-9.  Relati/e  center  ir radiance 
as  a  function  of  time  for  two  pulse 
forms  for  mild  blooming. 


Fig.  1-10.  Relative  center  irradiance 
as  a  function  of  time  for  two  pulse 
forms  for  strong  blooming. 


Fig.l-il.  Normalized  irradiance  ir 
fc.cal  plane  of  the  second  pulse  a  a 
function  :>f  radius  for  some  value'*  of 
energy  in  the  first  pulse. 


RAOUL  OtSTASCE  (cm) 


(2/  Propagation  of  Pulse  Trains 


The  effect  of  a  pulse  train  on  propagation  is  shown  in  Fig.  1-9  for  the  following  conditions: 
pulse  energy  -  80  kJ,  ranfc  -  5  km,  beam  quality  parameter  fi  2,  c  =  10-1>  cm-1,  and  a  hydro- 
dynamic  time  t^  =  140  psec  f>  r  a  single  pulse  20  psec  long  and  a  ..rain  of  five  pulses  with  the 
same  total  energy,  each  pulse  .8  psec  long.  The  pulse  train  shows  slightly  more  blooming  than 
the  single  pulse  as  a  consequence  of  the  hydrodynamic  results  shown  above. 

The  same  results  hold  for  a  pulse  with  strong  thermal  blooming,  shown  in  Fig.  1-10,  for  a 
pulse  energy  of  20  kJ  and  a  pulse  duration  Tp  =  100  psec  long.  Other  parameters  are  the  same 

as  in  the  previous  case.  ,  .. 

J.  Herrmai  n 

b.  A  Double- Pulse  Program 

A  propagation  program  has  been  written  which  handles  double  pulses.  The  first  pulse  (pre¬ 
cursor)  is  propagated  in  the  short-time  regime  (t-cube  aDpi  oximation).  The  energy  deposited 
during  this  pulse  is  used  to  heat  the  atmosphere  in  the  long-time  approximation.  If  the  second 
pul3e  follows  the  first  pulse  in  a  time  large  compared  to  the  hydrodynamic  time  it  will  encounter 
just  this  gas  density.  The  second  pulse  itself  again  is  in  the  short-time  regime.  No  dissipation 
mechanism  (wind,  sle  wing,  or  natural  convection)  is  included  between  the  two  pulses.  Infinite 
Gaussian  beams  are  vised. 

A  measure  of  the  strength  of  the  first  pulse  is  the  quantity 

E„  =  ae_aRk2E  (J/cm3)  , 

c  p 

where  E^  is  the  total  energy  in  the  pulse,  a  the  absorption  coefficient,  R  the  range,  and  k  the 
propagation  constant. 

In  Fig.  1-11,  the  irradiance  distribution  of  th-“  second  pulse  is  shown  as  a  function  of  radius 
for  some  values  of  the  energy  of  the  precursor  Ep.  The  first  pulse  for  these  cases  shows 
no  or  very  small  therm-  dooming.  The  cases  studied  are  for  the  following  conditions: 
range  R  =  3  km;  absorpu  coefficient  o  =  10  ^  cm  S  wavelength  =  2*/k  =  10.6  pm;  beam 
quality  parameter  /l  =  4;  transmitter  diameter  D  =  1  m,  infinite  Gaussian;  and  focal  area  A  = 

100  err.2.  /.  truncated  Gaussian  with  a  beam  quality  3  has  the  came  focal  area  and  therefore 
the  same  thermal  blooming. 

The  second  pulse  shows  large  thermal  blooming  if  the  energy  of  the  precursor  is  above  2  kJ, 

2 

or  (expr  essed  in  a  more  general  way)  if  Ec  is  above  3000  J/cm  . 

J  Herrmann 


c.  Four -Dimensional  Propagation  Code 

In  the  past  our  propagation  codes  have  been  limited  to  three  imensions:  the  CW  code  has 
three  spatial  dimensions,  and  the  pulse  codes,  assuming  circul  r  symmetry,  have  two  spatial 
dimensions  plus  time.  But  there  are  a  number  of  problems  which  develop  in  time,  but  without 
circular  symmetry.  For  instance,  the  propagation  of  a  CW  beam  through  a  stagnant  zo  ie  does 
not  lead  (under  our  hypotheses)  to  a  steady  state;  and  the  propagation  of  a  pulse  through  a  tur¬ 
bulent  medium  does  not  have  circular  symmetry.  Thus  it  is  clear  that  a  full  four-dimensional 
code  would  have  inte-esting  applications. 


j*. 
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FRACTION  OF  UNBLOOMEO  CENTRAL  INTENSITY 


Fig.  1-12.  Flux  within  a  circle  of  given  radius  in  the  focal  plane,  as  a  function 
of  radius,  for  four  initial  irradiance  distributions  with  and  without  blooming. 
Distributions  at  beginning  of  pulse,  before  blooming  starts.  Pbwer  -  8  x  10*>  W; 
pulse  length  -  1  nsec;  absorption  coefficient  -  1.5  x  10-3  cm'*;  range  -  4.0  jr,; 
X  =  15  pm;  and  neck  radius  -  3.127  cm  lor  all  cases.  Initial  conditions  for  four 
cases:  (a)  side  of  square  -  3.3  cm:  (b)  radius  of  uniform  circle  -  1 .94  crn; 
(c)  infinite  Gaussian,  l/e  power  radius  -  0.97  cm;  (d)  truncated  Gaussian,  l/e 
power  radius  -  1.37  cm;  radius  of  truncation  -  1.94  cm. 


Fig.  1-13.  Time  evolu  -on  of  a  modulated  pulse 
(long-time  approximation).  Times  in  units  of  r, 
defined  as  r  =  i.oc  T/2ok4(n  -  1)  Pe_orR. 
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Unfortunately,  a  simple  computation  also  makes  it  clear  that  a  straightforward  extension 
of  our  CW  code  into  the  time  dimension  by  use  of  a  standard  hydrodynamic  code  would  maka 
prohibitive  demands  on  both  computer  storage  ami  computing  time.  We  have,  therefore,  written 
a  code  that  extends  our  CW  code  by  a  single  time  step,  using  either  the  short-time  or  the  long¬ 
time  approximation  for  the  hydrodynamics.  Experience  with  our  circularly  symmetric  pulse 
codes  has  shown  that  these  approximations  are  valid  over  a  considerable  range:  the  short-timt 
approximation  gives  good  results  up  to  about  half  the  hydrodynamic  time,  while  the  long-time 
approximation  is  reasonably  good  after  2  to  3  hydrodynamic  times.  Thus,  it  should  be  possible 
to  cover  a  good  fraction  of  the  regimes  of  interest  for  pulse  propagation  using  the  new  cod  \ 

At  present,  this  code  does  not  include  the  capability  of  studying  the  stagnation  zone  problem, 
but  an  appropriate  modification  should  not  be  difficult.  The  computing  time  and  computer  stor¬ 
age  required  in  the  code  are  nearly  twice  those  for  the  corresponding  CW  computation. 

Using  this  code  we  have  compared  the  blooming  for  a  uniformly  illuminated  square  aperture 
t  1th  that  for  various  circularlv  symmetric  field  distributions.  The  distributions  were  chosen 
so  that  the  irradiance  in  the  focal  plane  with  nc  blooming  averaged  over  angles  was  nearly  the 
same  for  all  of  them  (Fig.  1-12).  The  conditions  correspond  to  a  laboratory  case.  When  blooming 
was  present  (short-time  approximation),  the  focal  plane  irradiance  was  again  nearly  the  same 
(Fig.  1-12),  showing  that  the  short-time  blooming  essentially  depends  only  on  the  mean  unbloomed 
behavior  in  the  immediate  vicinity  of  the  focus. 

L.  C.  Bradley 
J.  Herrmann 


d.  Brueckner  notabilities 

5  6 

Brueckner  and  Jorna  and  Kroli  and  Kelley  have  discussed  the  poBBibie  importance  of 
dynamic  instabilities  as  a  limitation  to  the  propagation  of  high-power  laser  beams.  These  arise 
from  the  coupling  of  the  laser  beam  wi*>.  entropy  or  sound  waves,  tnd  are  referred  to  as  stimu¬ 
lated  Rayleigh  and  Brillouin  scattering,  respectively.  These  analyses  were  linearized,  and  were 

7 

principally  concerned  with  infinite  plane  waves.  Rangnekar  has  made  some  numerical  studies 
of  instabilities  in  circularly  symmetric  collimated  beams,  using  a  ray-tracing  technique.  His 
results  confirmed  the  order  of  magnitude  of  Brueckner' s  analytic  formulas,  lr.  general  one  may 
say  that  for  CW  beams  these  effects  are  not  likely  to  be  important,  because  of  the-  comparatively 
small  irradiances  involved  and  because  of  the  smoe-ning  effect  of  convection.  But  for  pulsed 
beams  these  considerations  do  not  apply,  and  it  is  possible  to  get  sigr'ficant  amplification  of 
transverse  intensity  (or  phase)  fluctuation. 

Recently  indications  of  the  presence  of  this  effect  have  shown  up  in  some  of  our  nnioc  cal¬ 
culations.  The  observation  is  that  numerical  noise  at  the  Nyquist  frequency  appears  to  be  ampli¬ 
fied  in  nonlinear  propagation.  In  order  to  study  this  phenomenon,  I  have  used  our  propagation 
code  for  pulses  with  initial  conditions  corresponding  to  various  amplitude  modulations  of  the 
initial  beam.  An  example  is  shown  in  Fig.  1-13  where  the  initial  field  consisted  of  a  Gaussian 
plus  a  small  admixture  of  a  higher  order  eigenfunction  cf  circular  symmetry.  To  be  specific, 
the  field  was  taken  to  be  proportional  to 

[1  +  0.05  L12(r"  it*)  J  exp(-rZ/2a2) 

Here  Ln  is  the  Laguerre  polynomial  of  order  n.  It  can  be  seen  that,  as  time  increases,  the 
modulation  first  becomes  stronger,  and  then  largely  disappears.  In  the  present  case  the  overall 
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effect  of  the  modulation  on  the  blooming  is  small,  but  it  is  not  clear  that  this  will  be  true  for 
other  conditions.  In  particular,  it  is  possible  that  the  effects  of  turbulence  could  be  amplified 
by  this  mechanism.  Further  exploration  is  planned. 

L.  C.  Bradley 


e.  Propagation  of  Second  Pulse  Through  Fog  and  Aerosols 

We  can  consider  the  effect  of  pulse  beam  spread  due  to  a  fog  in  a  manner  analogous  to  our 
discussion  on  CW  beam  spread  due  to  aerosols  cited  previously  (cf  Sec.  1-A-l).  However,  for 
the  pulse  case  the  first  pulse  through  the  beam  is  greatly  attenuated  by  scattering  and  water 
vapor  absorption  until  the  absorbed  energy  exceeds  the  threshold  energy  necessary  for  vapori¬ 
zation  of  all  of  the  water 

E>Ethres  =  Ww<’d2l/4> 

where  hy  is  the  heat  of  vaporization,  Np  is  the  number  density  of  water  droplets,  D  is  the  beam 
diameter,  L  is  the  path  length  through  the  fog,  and  pw  is  the  density  of  liquid  water.  Assuming 
the  first  pu.se  has  sufficient  energy  to  clear  the  fog,  we  then  are  interested  in  how  well  the  next 
pulse  (following  by  a  few  microseconds,  for  example,  after  the  end  of  the  previous  pulse)  is  able 
to  propagate  through  the  cleared  channel.  The  extinction  coefficient  of  the  water  vapor  is  now 
considerably  less  than  that  of  the  water  drops,  and  we  can  make  a  preliminary  estimate  of  the 
beam  spread  due  to  multiple  scattering.  Because  after  a  few  microseconds  the  index  spheres 
have  expanded  to  a  size  much  larger  than  the  wavelength  (cf.  Sec.  1-03),  we  may  use  the  geo¬ 
metric  optics  approximation  for  multiple  scattering  through  the  ensemble  of  index  spheres.  As 
shown  before  for  the  CW  case,  the  mean  squared  deflection  is 

yrZ  =  4/9  L3Nr*nZ(t)  (An/n)2  . 

For  the  pulsed  case,  we  assume  that  the  water  droplet  heats  up  to  the  critical  temperature, 

T  =  374"K,  and  heat  is  then  transferred  rapidly  to  the  air  into  a  rapidly  expanding  heated  air- 
vapor  sphere.  The  composition  of  the  vapor  sphere  <s  predominantly  heated  air  whose  velocity 
we  hav»  modelled  empirically.  We  assume  that  the  change  in  index  of  refraction  is  uniform 
throughout  this  index  sphere  and  is  then  given  by 

3 

r 

An  =  0AT 

aJ(t) 
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An  appropriate  distribution  for  fog  radii  is  given  by 


dN_  =  N 
r  o 


<rm> 


dr 


where  r  is  the  peak  of  the  distribution,  so 
m 


^total 


'Ll 


1/2  ,  (7!)1/2N1//Zi1'/2r^ATL3'/2 


a2(t)  n 


If  we  compare  this  with  the  corresponding  diffraction-limited  spreading  for  the  pulsed  case  in 
fog,  we  have 


"total 

^diff  fog 

pulsed  beam 


(7i)1/2ii/2L1,/2N  *^2r  30ATD 


a  (t)  nX 


For  values  of  cpy  =  2  j/g'K,  cpa  =  1  J/g*K,  py  =  1  s/cm,9  Tq  =  300*K.  ag  =  0.25  cm3/g. 

L  =  1  km,  D  =  2  cm,  rm  =  0.6  pm  (e.g.,  as  measured  by  Harris9),  and  n  =  1,  w^  "Kow  in 
Fig.  1-14  the  development  in  time  of  the  rela.ive  beam  spread  for  three  different  fog  droplet 


Fig.  1-14.  The  rms  total  defle  "tion  normalized  to  tiv*  diffraction  spread 
for  a  pitlsed  10.6-pm  laser  beam  of  uniform  intensity  and  2-cm  diameter 
through  a  fog  (with  a  peak  in  particle  radius  distribution  of  rm  =  0.6  pm) 
of  a  path  length  of  1  km  for  three  water  droplet  densities. 


densities.  A  caution  here  is  that  fogs  can  range  to  distributions  with  a  value  of  r^  as  high  as 
5  pm  and  since  the  deflectio..  scales  as  r^,  the  deflections  can  become  considerably  more 
severe  than  those  indicated  by  Fig.  1-14.  The  empirical  value  of  a(t)  described  in  Sec.  I-C-3 
of  this  report  has  been  used. 
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TABLE  1-1 

COMPARISON  OF  PHYSICAL  PARAMETERS 
OF  CHIEF  AEROSOL  COMPONENTS 


Boiling  Point 
(°C) 

Specific  Gravity 

C 

P 

(J/gm  °K) 

c 

p  p 

aei  aer 

c  p 

p  .  w 

water 

AT 

aer 

AT 

w 

Carbon 

4827 

2.25 

0.17(0.3) 

1.9 

S:°2 

2230 

2.2 

0.22 

1.54 

NaCI 

1413 

2.16', 

0.21 

0.91 

If  we  assume  that  in  the  pulse  regime  a  heated  aerosol  acts  very  much  like  a  water  droplet, 
evaporating  rapidly  with  a  corresponding  rapidly  expanding  sphere  with  approximately  the  same 
time  evolution,  we  have  for  the  relative  rms  beam  spread 


ytotal 


ydiff 


aerosol 
pulsed  beam 


r1/?,L1/<JN1/2(r  r  .  )3/2/i ATD 
o  max  min 

- - - 

a  (t)  nX 


-4  3 

where  we  have  assumed  again  for  aerosols  the  distribution  1  =  3r  N  r  ,  dr  between  radii  of 

r  o  min 

rmin  anc*  rmax'  ’ms  deflecti°n  due  to  aerosols  relative  to  that  due  to  fags  for  the  same 

concentrations  is  then  given  by 


,  >1/23  AT  c  p 

yaerosol  _  f  rrrtaxrmin^  j  1  “er  paer  aer 

yfog  l  rm  J  (7'.)3^2  ^^wa^pv*1  w 


In  Table  I-l  we  list  three  of  the  dominant  components  of  aerosols  with  their  boiling  points  (hence 
approximate  critical  temperature  change)  and  other  relevant  physical  constants  and  compare  tne 
factor 


c  p  AT 
paer  aer  aer 

c  a  AT  . 
p\r  w  water 

Multiplying  the  last  column  by  0,02  would  give  approximately  the  spread  due  to  aerosols  relative 
to  a  fog  of  the  same  concentration  of  scatterers.  These  results  give  us  an  initial  indication  of 
the  significance  of  beam  spread  due  to  fogs  and  aerosols. 


f.  Propagation-Through-Fog  Code 

In  order  to  solve  a  more  accurate  formulation  of  the  beam  spread  of  the  second  pulse  through 
a  cleared  log,  we  intend  to  modify  the  Bradley-Herrmann  propagation  code30  to  include  the  scat¬ 
tering  effects  due  to  the  vapor  spheres  formed  around  the  water  droplets  (r  -  he  index  spheres 
that  form  around  Irradiated  aerosol  particles). 

In  the  Bradley-Herrmann  code,  the  effect  of  the  medium  is  simulated  by  a  phase  change  of 
the  field  propagated  through  a-small  axial  increment.  To  employ  the  code  we  need  to  find  if  for 
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sufficiently  low  spatial  frequencies  the  sum  of  scattered  and  incident  field  can  be  described  ac- 
curatelv  by  a  simple  phase  change,  i.e.,  at  each  point  in  a  plane  of  observation  beyond  the  scat- 
tereis  and  perpendicular  to  the  propagation  direction  of  the  incident  beam, 


E.  +  E 
inc  scatt 


E.  e 
inc 


ik  / n(z)dz 


where  kQ  is  the  optical  wave  number  and  n(z)  is  the  index  of  refraction.  For  index  changer  not 
very  large  the  exponential  can  be  expanded  to  first  oroer  and  the  incident  field  terms  subtracted 


E  ..  =  F.  S(0) 
scatt  inc 


exp[(-R  +  z)  ik  J 


ik  R 
o 


a  E 


.  ik  r 

inc  o  J 


n\z)  dz 


We  then  wish  to  construct  the  given  integrated  index  of  refraction  (or  its  corresponding  spatial 
Fourier  transform  in  the  plane  of  observation)  and  check  the  validity  of  our  low  spatial  frequency 
approximation. 

The  Fourier  transform  r(k)  of  the  integrated  phase  change  through  the  sphere  is  defined  by 

k  f  ndz  =  \  v(k)  e1^  rdZk 
o  J  J 


We  wish  to  compare  this  with  the  Fouiier  transform  of  the  scattered  wave  v'(k)  defined  by 


S,(0)  exp  [(-R  +  z)  ik  I  p  jt  •• 

-  J- - 0 - ^  ■  j  .■(»  .'k  r. 


d2k 


The  geometry  of  this  equation  is  shown  in  Fig.  1-15  where  R  is  the  vector  from  the  center  of  the 
scatterer  to  a  point  of  observation  in  the  plane  perpendicular  to  the  incident  direction,  z  is  the 
vector  along  the  in  ident  direction  from  the  scatterer  to  the  incident  plane,  kQ  id  the  optical  wave 


Fig  1-15.  View  of  the  geometry  of  scattered 
radiation  problem  for  a  single  fog  droplet. 


number,  and  Sj(©)  are  the  scattering  amplitude  functions  for  the  perpendicular  and  parallel  elec¬ 
tric  field  components  as  a  function  of  the  angle  &  between  z  and  It.  Also,  k  is  the  spatial  wave 
vector  in  the  plane  of  observation  and  T  is  the  vector  from  the  intercept  of  z  and  the  plane  to  the 
point  of  observation  in  the  plane.  The  spatial  transform  for  a  single  particle  is  then  given  by 


•'’(k) 


dr  r  JQ(kr) 


cos©  S^e)  exp ((-  e0Z--§  +  z)  ik,! 

- " - n. - 

o 


(l~i) 


IS 


For  times  greater  than  a  few  microseconds,  the  vapor  spheres  become  much  larger  than  the 
wavelength  of  the  incident  lip1-*.,  and  the  total  scattering  component  from  a  single  sphere  is 
much  less  than  the  incident  .eld  on  the  sphere,  so  it  is  appropriate  to  employ  the  Rayleigh-r,ans 
approximation  for  the  scattered  field*  * 

S.<0)  =  ik  a  An  - - vjy  Jw,(2k  a  sin  1/2  6) 

}  °  (2kQa  sinl/20)3/Z  3/2  °  |cos0  .  (1-2) 

In  Eq.  (1-2)  the  respective  upper  and  lower  scattering  amplitudes  are  for  electric  field  components 
perpendicular  and  parallel  to  the  reference  plane,  and  An  is  the  uniform  index  change  within  a 
sphere  of  radius  a  which  is  expanding  in  time.  By  using  trigonometric  identif.es  ar.d  variable 
substitutions,  we  can  rewrite  (1-1)  as 

*  ,i/4*°' 3/2,172'.  X  [I  ds  kz  8)  J3/2<,/2  koaa) 

to  (K  cl)  *  1  •  O 

o 

Xs  1//2  exp[-(ikQZB2  +bs2)ij  (1-3) 

where  we  have  assumed  that  s  «  1  (the  approximations  involved  are  within  SO  percent  accuracy 
even  for  0  =  »/ 3  or  s  »  1)  anf  introduced  the  expression  exp(-bs2)  to  limit  contributions  for 
’.arge  s  (i.e.,  b  must  be  at  le  st  comparable  to  1).  The  Introduction  of  the  exponential  limits 
contributions  from  the  integral  for  large  s  which  only  occur  as  a  result  of  our  approximations, 
thus  allowing  the  integral  to  converge.  Additionally,  only  the  small-angle  scatter  contributions 
interest  us  since  large-v  gle  scatter  produces  an  effective  attenuation.  The  integral  (1-3)  can 
be  done  exactly  and  the  result  is: 


k2a3An 
o _ 


«  2  2 

Tk^zTB  l  [if&zVb)]  F(_m*  "m;  5/2:  77)  (1‘4) 


where  F  is  the  hypergeometric  function. 

For  large  values  of  z  and  for  values  of  a  not  too  large,  both  u  and  ou  will  be  much  less  than  1, 
and  we  can  curtail  the  expansion  to  the  leading  terms  in  u  and  retain  terms  only  to  second  order 
in  a.  We  then  need  for  our  analysis  the  spatial  spectrum  <V(k)  v'*  (k)>,  which  is  averaged  over 
z  since  the  axial  position  of  the  particle  is  a  random  variable.  Retaining  terms  only  to  fourth 
order  in  k  are  find 


<s'(k)  s'*(k)> 


"  ,7 


H 


If  w.  compare  this  to  the  corresponding  geometric  optics  solution  for  light  passing  through  a 
uniform  sphere ,  we  find  the  spatial  spectrum  to  be 


<*(k)  a*  (k)> 


2  6  * 

km  An" 
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The  agreement  to  first  order  allows  us  to  use  (1-5)  in  place  of  (1-4)  in  our  propagation  code. 

We  next  must  incorporate  the  single-particle  spectrum  into  the  spectrum  of  a  density  of 
randomly  distributed  such  particles,  and  we  shall  then  be  able  to  begin  our  computer  Simula  ion 
of  the  distortion  of  the  second  pulse  through  an  evaporated  fog  channel.  Our  initial  assumptions 
for  the  multiparticle  system  are  that  the  position  of  the  particles  are  uncorrelated,  which  yields 
the  almost  white  noise  spectrum.  However,  we  anticipate  that  we  shall  eventually  require  cor¬ 
relation  of  particle  position  to  account,  for  example,  for  tufts  of  fog. 

II.  Granek 
I..C.  Bradley 


f.  Propagation  of  a  Focused  Beam  Through  an  Absorbing  Medium 


This  section  considers  ti.e  extent  to  which  the  atmosphere  CO^  can  be  described  as  a  satu¬ 
rable  medium  lo  10  6 -pm  radiation  with  a  phenomenologica’  absorption  coefficient  given  by 

«  =  oo/(l  +  i/y  . 

For  a  simple  two-level  system  of  the  type  described  by  Selden12  (Fig.  1-16),  the  interacticn  with 
radiation  will  cause  a  transient  change  in  level  populations  followed  by  a  steady-state  absorption. 
The  time-dependent  coupled  equations  describing  such  a  simple  system  are 


_dl_ 

dZ 


-ol(n2  -  rij) 


hr  <n2  nl'  n2 /T 


o 


dn 


VI  I  ‘l  ■» 

~dT  =  +  hr(nZ  ~  V  +  n2^To 


Here,  1  =  intensity,  n^  =  density  of  the  upper  level,  n^  =  density  of  the  lower  level,  a  =  absorption 
cross  se  on,  and  Tq  is  the  relaxation  rate  of  the  upper  level. 

|il-i-wtl| 


Fig.  1-16.  Illustrated  is  a  simple 
♦wo-lev.-l  system  which  absorbs 
laser  radiation. 


LASER  TRANSITION  LINE 
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After  the  initial  transient  radiation  edge,  a  steady  statr  will  evolve  which  is  governed  by  the 
incident  intensity  and  the  effectiveness  of  the  relaxation  of  the  upper  level  with  relaxation  time 
Tq.  In  such  a  steady  state  me  temporal  rates  of  change  of  the  populations  are  zero.  By  sub¬ 
stituting  the  second  rate  equation  into  the  first,  the  rate  of  change  of  the  intensity  is  given  in 
terms  of  the  upper  level  only.  Now  using  the  second  equation  and  the  condition  that  the  total 
population  of  both  levels  is  some  constant  value  n,  the  following  equation  results: 

_dl  -on! 

dz  (i  +  i/y 

where 


I  = 


hv 

2ffTo 
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Fig.  1-17.  Illustrated  is  a  three-level 
system  which  approximates  absorption 
in  the  atmosphere  by  CC^. 


T 


1 


Fig.  1-19.  A  plot  of  P(Z1/P(0)  versus  l(0)/l.  with  A.  *  104,  A  .  =  10z, 
f  »  105  cm,  and  o  *  lO’*  8  “  out 
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Considering  a  three-level  system,  illustrated  in  Fig.  1-1 7,  the  same  procedure  will  result 
in  the  following  exp-,  ession: 


dZ  "  2(1  +  1/1  ) 


where  now 


i  2  hr 

*8  =  (3T2  +  Tj)  a 

Here  Tj  is  the  relaxation  rate  for  the  second  level  to  the  first,  and  is  the  relaxation  rate 
of  the  third  level  to  t  ie  first  by  collisions. 

For  the  propagation  of  10.6-pm  radiation  through  the  atmosphere,  this  simple  three-level 
system  will  be  used  estimate  the  saturation  intensity?  The  upper  laser  level  in  CO^  is  quickly 
relaxed  by  which  remains  excited  for  the  duration  of  the  pulse.  The  rate  at  which  the  lower 
laser  level  is  populated  depends  on  T-V  collisional  excitation  of  the  lowest  vibrational  level. 
Water  vapor  contributes  heavily  to  this  process  and  Tj  is  given  by 


where  T^  is  due  to  other  molecules. 

For  a  focused  beam,  the  usual  phenomenological  expression  for  the  rate  of  change  of  inten¬ 
sity  with  distance  must  be  modified  to  account  for  the  increase  in  the  intensity  due  to  a  fixed 
power  being  concentrated  over  a  progressively  smaller  area.  This  gives  an  additional  term  as 
shown  below 

dl  8A  i  -al 

dZ  6Z  A  '  (1  +  I/lg) 

By  redefining  I  =  yIsAo/A(Z)  and  letting  Za  =  p,  the  equation  achieves  the  compact  form 


This  yields  the  following  formal  solution 


y  =  y(0)  exp 


° 


Thus,  the  value  of  y  at  p  is  related  at  p  +  6p  by  the  approximate  relation 


y(p  +  6p )  “  yip)  exp 


'U  +  yip)  Anl 

A(p) 


By  using  this  simple  iteration,  a  numerical  evaluation  of  P(Z)/P(0)  =  yfaZ)/y(0)  is  extremely 
efficient  resulting  in  valuet  of  the  type  shown  in  Figs.  1-18  to  1-20.  These  figures  are  represen¬ 
tative  of  a  laboratory  situation  (a  10-cm  beam  focused  to  a  3-mm  spot  in  5  meters)  and  two  field 
situations  (a  1 -meter-diameter  beam  focused  to  a  10-cm  spot  over  1  km  and  a  collimated  10-cm- 
diamoter  beam) .  \ 


D.  L.  Moone^ 
H.  Kleiman 


J.  E.  Lowder 
R.  W.  O'Neil 
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Fig.  1-20.  The  graph  of  ?(z)/P(0)  with  constant  area  of  102  cm2  and  q 
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Fig.  1-21.  Geometry  for  saturation  calculations. 
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Fig-  1-22.  Plot  of  ratio  of  the  incident  power  to  the  transmitted  power  as  a  function 
of  the  ratio  of  the  transmitted  intensity  to  the  saturation  intensity. 
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2.  Experiments 


a.  Pulsed  Thermal  Blooming  and  Saturation  of  CO ^  at  10.6  pm 

The  experiments  described  in  the  previous  Optics  Research  Report'1  have  been  analyzed. 

The  results  for  absorption  by  a  broadband  absorber  (propane)  causing  thermal  blooming  within 
a  hydrodynamic  time  have  been  submitted  for  publication.11 

A  number  of  experiments  were  also  done  in  which  the  primary  absorber  was  500  torrs  CO^ 
in  250  torrs  argon.  The  primary  physical  process  we  were  inter. sied  in  examining  here  was 
the  saturation  of  the  CO^.  For  pulse  lengths  of  the  type  used  (~i0  psec)  it  is  justified  to  use 
steady-state  kinetics  in  describing  the  saturation.  Hence,  we  describe  the  absorption  coefficient 
phenomenologically  by: 

a  =  a0/l  +  (i/lg)  . 


Here,  aQ  is  the  linear  absorption  coefficient. 
Then  using  the  geometry  shown  in  Fig.  1-21 


S-- 


”op 


1  + 


?  d2<Z>  lS 


where 


Here,  P  -  power,  A.  =  wavelength,  Z^  the  focal  length,  Z  the  distance  as  measured  from  the  focal 
point,  and  ig  the  saturation  intensity. 

One  can  solve  this  differential  equation  numerically.  If  one  solves  for  P(Z),  then  P(Z  =  0)/ 

P(Z  =  Z,)  is  the  fraction  of  the  power  transmitted  through  the  absorbing  cell.  Using  1  as  an 
I  s 

adjustable  parameter  we  have  fit  curves  of  this  type  to  the  data  as  shown  in  Fig.  1-22.  For  the 
conditions  of  our  experiment  (a  =  1.8  X  10  ^cm  *,  Z^  =  490  cm),  we  obtain  1  s  1  to  2  X  10^W/cm2. 

We  are  currently  working  on  more  sophisticated  theory  detailed  in  Sec.  U-A-4  to  explain  the 
magnitude  of  the  saturation  parameter. 

H.  Kieiman 
R.  W.  O'Neil 
J.  E.  Lowder 


b.  Laser-induced  Air  Breakdown  for  1.06-pm  Radiation 

Experiments  have  been  done  at  1.06  pm  for  the  breakdown  threshold  in  clean  air  and  in  car¬ 
bon  aerosol  as  a  function  of  spot  size.  A  Q-switched  neodymium  glass  iaser  providing  6  J  in  a 
100-nsec  pulse  was  used  for  these  experiments. 

The  breakdown  threshold  (50  percent  probability  of  breakdown)  in  clean  air  was  measured  in 
a  chamber  with  highly  filtered  air  for  four  spot  sizas.  The  spot  si’e  of  10  2  cm  was  the  largest 
value  for  which  clean  air  breakdown  could  be  obtained  (on  *10  percent  of  the  shots)  with  the  avail¬ 
able  pulse  power.  The  difference  between  the  laboratory  air  threshold  and  the  clean  air  value 
was  observable  only  for  the  largest  spot  where  the  laboratory  air  threshold  was  *0.5  thai  in 
clean  air.  The  data  are  shown  in  Fig.  i-23. 
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SPOr  SIZE  It"’  diomettr  cm) 

Fig.  1-23.  Breakdown  threshold  in  clean  air  as  a  function  of  spot  size.  Solid  curves 
are  from  Ref.  17  with  A  =  D/7  a;id  Tp  =  100  nsec.  The  error  bars  indicate  threshold 
uncertainty.  The  relative  accuracy  i3  «<A20  percent. 


PARTICLE  DIAMETER  turn) 

Fig.  1-24.  Threshold  in  carbon  aerosol  as  a  function  of  particle  size. 
The  value  at  50  pm  was  measured  for  single  particles  in  a  270-pm 
spot.  The  values  for  smaller  particles  were  obtained  for  a  100 -pm 
spot  in  a  carbon  dust  atmosphere  where  the  particle  size  indicates 
the  estimated  value  for  the  largest  particle  in  the  focal  volume.  The 
relative  accuracies  for  the  thresholds  are  *20  percent. 
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The  results  for  clean  air  indicate  somewhat  higher  thresholds  than  have  been  reported  pre- 
14-16 

viously.  There  are  several  factors  which  can  account  for  this.  For  small  spot  sizes  the 

threshold  js  verv  sensitive  to  diffusion  length  which  depends  on  the  intensity  distribution  in  the 

focused  spot,  while  for  larger  spots  the  threshold  can  be  lowered  by  the  presence  of  naturally 

1 7 

occurring  laboratory  dust  as  will  be  discussed  below. 

The  solid  curves  in  Fig.  1-23  are  theoretical  values  obtained  from  Ref.  18  for  =  100  r.sec 
and  a  diffusion  length  A  =  D/7.  The  different  curv.  „  ''orrespond  to  the  inclusion  in  their  calcula¬ 
tions  of  various  photoelectric  processes  requiring  n  or  fewer  photons,  n  =  0  corresponds  to  the 
classical  microwave  theory. 

These  results  suggest  that  one  and  possibly  two  photon  processes  are  involved  in  1.06-pm 
breakdown  but  the  dominant  mechanism  appears  to  be  cascade  ionization.  The  k  ^  scaling  of 
microwave  theory,  which  is  valid  at  10.6  pm  (Ref.  17),  is  just  becoming  invalid  at  1.06  pm.  A 
similar  conclusion  was  arrived  at  .com  experiments  involving  1.06  gin  and  shorter  wave¬ 
lengths.16'19 

The  effects  of  particles  on  the  threshold  were  investigated  using  the  techniques  described  in 
Refs.  17  a-d  20.  Single,  50-pm  carbon  particles  were  positioned  at  the  center  of  a  270 -pm  fc- 

o  ,  2 

cused  spot.  The  threshold  for  this  case  was  »5  x  10  W/cm  or  »0.02  of  the  clean  air  threshold. 
The  threshold  was  then  measured  in  a  carbon  aerosol  for  a  100-pm  spot.  The  measured  values 

The 

dust  content  in  the  beam  can  be  characterized  by  the  size  of  the  largest  particle  in  a  focal  vol¬ 
ume.  The  threshold  data  are  shown  in  Fig.  1-24  for  the  corresponding  particle  size.  The  dashed 
curve  represents  10*'  times  the  similar  threshold  measured  at  10.6  pm  (Ref.  17).  The  thrf  sholds 
meast.red  for  the  100-pm  spot  in  a  dust  environment  can  be  compared  directly  with  the  corre¬ 
sponding  clean  air  threshold  shewn  in  Fig.  1-23:  these  being  from  0.1  to  0.4  times  the  clean  air 
value.  It  should  be  noted  that  the  values  in  the  dusty  air  are  comparable  to  the  previously  re¬ 
ported  thresholds  in  air. 14  1  J 

The  primary  conclusions  from  this  paper  are  that  (1)  cascade  ionization  is  the  dominant 

breakdown  mechanism  at  1 .06  pm  with  contributions  from  one-  and  two-pl.oton  processes, 

-2 

(2)  threshold  measurements  in  air  for  spot  sizes  larger  than  “10  cm  can  be  lowered  by  the 
presence  of  laboratory  dust,  and  (3)  the  fractional  reduction  is  comparable  to  those  observed  at 
10.6  pm  p.  g  Lencioni 


for  light  and  heavy  dust  conditions  were,  respectively,  »2.5  x  lo11  and  6  x  lo10  W/cm2. 


c.  Dynamics  of  Particle-Induced  Air  Breakdown  in  A  Laser  Beam 

Further  experiments'  on  air  breakdown  triggered  by  single  50-pm  carbon  particles  were 

done  to  study  the  growth  dynamics  of  the  breakdown  within  the  laser  beam  as  a  function  of  laser 
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intensity,  and  the  absorption  length  in  the  breakdowr  plasma. 

The  experimental  apparatus  is  shown  in  Fig.  1-25.  The  laser  beam  was  provided  by  a  5  J 
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per  pulse  Febetron-CC^  laser  operated  in  the  short-pulse  mode.  An  image  converter  camera, 
with  exposure  times  down  to  10  nsec,  was  used  tc  obtain  time -resolved  information  on  the  growth 
of  the  breakdown  in  the  beam.  The  transmitted  beam  was  reimaged  and  magnified  9:1  at  a  pair 
of  photon  drag  detectors  to  obtain  the  radial  and  temporal  dependence  of  the  transmitted  beam. 

Figure  1-26  shows  a  series  of  image  converter  frames  which  illustrate  the  growth  of  a 
breakdown  at  near  threshold  intensity.  The  dashed  white  lines  mark  the  e  1  diameter  of  the 
beam  which  was  1  mm.  The  laser  beam  was  incident  from  the  left.  The  four  frames  were  taken 
on  four  separate  shots  at  the  indicated  exposure  times  relative  to  the  start  of  the  pulse.  The 
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Fig.  1-25.  Experimental  apparatus  to  study  breakdown  dynamics.  Time-resolved 
information  on  the  growtn  dynamics  was  obtained  with  the  image  converter  camera. 
Particles  were  placed  at  a  spot  in  the  focal  volume  with  the  activated  particle  dropper 
shown  in  the  insert. 
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Fig.  1-26.  Series  of  image  converter  frames  illustrating  breakdown 
growth  at  near-threshold  intensity. 
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Fig.  1-27.  Breakdown  dynamics  at  five  times  threshold  intensity. 
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breakdown  initiates  at  near  the  peak  laser  intensity  and  grows  radially  and  axially  as  a  volume 
absorbing  plasma  The  axial  velocity  is  slightly  larger  than  the  radial  velocity  which  had  an 
initial  value  »1.4  x  10  cm/sec.  The  velocities  were  observed  to  decrease  by  »25  percent  after 
«300  psec  due  probably  to  the  natural  decay  in  laser  intensity.  Note  that  the  particle  remains 
visible  in  each  of  the  frames. 

At  higher  laser  intensities  the  character  of  the  breakdown  growth  changed  abruptly.  This 

is  illustrated  in  Fig.  1-27  where  the  intensity  was  increased  above  threshold  by  a  factor  of  5.  In 

these  frames  a  white  dot  has  been  inserted  to  indicate  the  position  of  the  particle.  Note  that  for 

these  higher  intensities  the  breakdown  develops  as  a  thin  absorbing  front  which  propagates  back 

towards  the  laser  and  grows  radially.  Note  also  the  increase  in  velocities.  This  type  of  break- 

23 

down  growth  is  characteristic  of  the  LSD  described  by  Raizer. 

The  radial  and  axial  velocities  obtained  from  data  of  this  type  are  shown  in  Fig.  1-28  as  a 
function  of  intensity.  At  lower  intensities  the  velocities  increase  rapidly  with  intensity  roughly 
proportional  to  1.  At  higher  intensities  the  axial  velocity  varies  roughly  a«  1*  ,  which  is  in 

agreement  with  planar  LSD  wave  theory. 
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Fig.  1-28.  Axial  and  radial  velocities 
of  breakdown  plasma. 


Fig.  1-29.  Absorption  depth 
in  breakdown  plasma. 


The  absorption  of  laser  intensity  within  the  breakdown  for  a  spot  diameter  of  1  mm  was 
measured  directly  by  reimaging  the  focal  spot  with  magnification  on  two  photon  drag  aetectors: 
one  giving  the  power  transmitted  through  the  central  disk  of  radius  r  =  0.028  cm  and  the  other 
giving  the  power  through  a  concentric  annulus  with  inner  radius  r^  and  outer  radius  0.067  cm. 
The  absorption  depth  within  the  breakdown  was  estimated  by  comparing  the  luminous  thickness 
with  the  transmitted  intensity.  The  res\  Us  are  shown  in  Fig.  -29. 


D.  E.  Lencioni 
L.  C.  Pettingiil 
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1.  Double -Pulse  Pressure  and  Impulse  Measurements 


Measurements  of  surface  pressure  and  impulse  delivered  to  aluminum  and  carbon  targets 
are  reported  in  which  two  laser  pulses  spaced  by  a  variable  time  delay  were  used.  The  meas¬ 
urements  are  presented  in  a  normalized  form  to  show  the  effect  of  the  first  pulse  on  the  peak 
pressure  and  impulse  delivered  by  the  second  pulse. 

The  experimental  arrangement  is  shown  in  Fig.  1-30.  The  two  CO^  e-ueam  iasers  are  de¬ 
scribed  in  Refs.  24  and  22  and  were  mounted  in  a  lead-lined  room  as  indicated  in  Fig.  1-30.  The 
3-cm-diameter  laser  beams  were  made  parallel  wi;h  about  a  i-cm  separation  between  them  and 
propagated  through  the  optical  system  shown  in  Fig.  1-30.  The  two  beams  were  focused  onto  the 
same  spot  at  the  target  surface  with  either  a  15-cm  or  a  28-om  focal  length  mirror.  The  effect 

of  the  convergence  angle  of  the  two  beams  is  shown  in  an  open-shutter  photograph  of  the  two 

24 

laser-produced  plasmas  in  Fig.  1-31.  One  laser  was  operated  in  the  long-pulse  mode  with  an 

energy  of  about  10  J  and  a  pulse  length  which  varied  from  40  to  50  psec.  The  second  laser  was 

22 

operated  in  the  TEA-pulse  mode  with  an  energy  of  about  7  J  and  a  pulse  length  of  from  2  to 
3  psec.  Typical  pulse  shapes  for  these  two  lasers  are  shewn  in  Fig.  1-32.  The  focal  spot  sizes 


Fig.  1-32.  Typical  pulse  shapes 
of  (a)  short-pulse  laser,  (b)  long- 
pulse  laser  measured  with  photon 
drag  detectors. 
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Fig.  1—3 3.  Peak  pressure  on  an  aluminum  surface  due  to  irradiation 
by  short  pulse  after  long-pulst  precursor  normalized  by  the  peak 
pressure  due  to  the  short  pulse  in  the  abcence  of  a  precursor.  The 
open  symbols  indicate  second-oulse  peak  power  densities  of  about 
7  x  10»  W/cm2  and  the  solid  symbols  indicate  power  densities  of 
about  2  X  to**  W/ci.i2. 
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Fig.  I -34.  Impulse  delivered  to 
0.38-cm-diameter  aluminum  and 
carbon  targets  as  a  function  of 
incident  single-pulse  energy.  The 
error  bars  denote  the  estimated 
maximum  error  of  the  impulse 
neaSurement. 
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used  in  obtaining  peak  power  densities  were  obtained  by  attenuating  the  incident  beam  with  fine 

perforated  screens  ard  observing  the  variation  of  the  resulting  burn  areas  on  thermally  sensitive 
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paper  Peak  power  densities  incident  on  the  target  were  a  few  x  10  W/cm  foi  the  long-puise 

8  2 

laser  and  a  few  x  10  VV/cm  for  the  short-puise  laser. 

Relative  pea1'  pressure  measurements  in  the  focal  spot  we'e  made  with  a  PCB/Piezotronics 
quartz  pressure  gauge  which  was  covered  with  a  1-mil  sheet  of  pure  aluminum.  The  thin  alumi¬ 
num  sheet  was  attached  to  the  0. 38-cin-diameter  diaphragm  of  the  quartz  gauge  with  a  thin  layer 
of  silicone  grease.  The  peak  pressure  was  measured  first  for  the  short  pulse  alone,  then  the 
pressure  resulting  from  the  short  pulse  preceded  by  the  long  puise  was  measured  for  various 
delay  times.  The  measured  single-pulse  peak  pressures  ranged  from  about  10  bar  for  the  long- 
pulse  laser  10  about  40  bar  tor  the  short-pulse  laser.  These  measurements  are  lower  limits  to 
the  actual  peak  pressures  at  the  target  surface  due  to  the  i-psec  response  time  of  the  pressure 
gauges,  but  the  measurements  do  provide  a  gsod  relative  measure  of  the  actual  single-  and 
double-pulse  target  pressures.  The  normalized  peak  pressure  measurements  for  the  short  pulse 
preceded  by  the  long  pulse  are  shown  in  Fig.  1-31  as  a  function  of  delay  time.  These  results  3how 
the  peak  pressure  caused  by  the  short-pulse  laser  is  reduced  to  about  10  percent  of  its  normal 
value  right  after  the  long-pulse  precursor,  and  it  takes  a  delay  of  the  firing  time  of  the  short 
pulse  of  f-om  1  to  2  msec  for  the  peak  pressure  to  recover  to  its  single-pulse  value. 

The  impulse  delivered  to  0.38-cm-diameter  aluminum  and  carbon  targets  was  measured  with 
a  sensitive  ballistic  pendulum.  The  impulse  was  found  from  the  mass  of  the  pendulum  and  the 
maximum  amplitude  of  its  swing  which  was  measured  with  a  photodiode-light  source  assembly. 
The  pendulum  was  calibrated  botn  by  measuring  the  period  to  determine  its  effective  length  and 
by  imp’’ "ting  it  with  small  bodies  of  known  momentum.  The  single-pulse  impulse  measurements 
for  both  the  long-  and  short- pulse  lasers  are  shown  in  Fig.  1-34.  The  impulse  delivered  by  the 
short-pulse  laser  was  found  to  vary  roughly  as  the  cube  root  of  the  incident  energy  with  a  little 
more  impulse  being  delivered  to  carbon  targets.  For  a  fixed  pulse  shape  and  focal  spot  size, 
the  peak  power  density  incident  on  the  target  is  proportional  to  the  energy  in  the  pulse.  Thus, 
the  energy  dependence  shewn  in  Fig.  1-34  is  in  agreement  with  results  presented  in  Refs.  26,  27, 
and  28  which  show  the  impulse  delivered  to  small  targets  is  proportional  to 

1  a  A^giq)1/3  , 

where  AT  is  the  target  area,  Dg  is  the  focal  spot  diameter,  and  q  the  incident  power  density. 
The  impulse  delivered  by  the  long-pulse  laser  was  found  to  be  much  higher  for  carbon  them  for 
aluminum  targets  as  is  shown  in  Fig.  1-34.  The  impulse  delivered  by  the  long-pulse  laser  was 

strongly  influenced  by  blowoff,  and  since  carbon  has  a  much  higher  surface  absorptivity 

29  29 

(o  “  85  percent)  for  10.6-pm  radiation  than  aluminum  (a  =>  3  percent),  the  impulse  was  cor¬ 
respondingly  greater.  The  relationship  between  impulse  find  energy  for  each  set  of  data  shown 
in  Fig.  1-34  was  used  to  determine  the  impulse  delivered  for  the  case  of  each  pulse  alone.  The 
lasers  were  then  fired  in  the  double-pulse  mode,  with  the  long-pulse  laser  as  the  precursor, 
and  the  total  impulse  transferred  to  aluminum  and  carbon  targets  was  measured  as  a  function  of 
delay  time.  The  results  are  presented  ir.  Fig.  1-35  where  the  total  impulse  measured  in  the 
double-pulse  mode  is  norma  -ed  by  the  sum  of  the  impulses  which  would  have  been  delivered  if 
each  laser  had  been  fired  separately.  The  results  in  Fig.  1-35  show  the  total  impulse  deli/ered 
to  the  targets  varies  from  50  to  80  percent  of  the  value  to  be  expected  for  single  firings.  The 
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Fig.  1-35.  Total  impulse  delivered  to  0.38-cm-diameter  carbon  and  aluminum 
targets  for  double-pulse  firing  normalized  by  the  impulse  delivered  by  each 
pulse  when  fired  separately.  These  data  were  taken  with  the  long  pulse  as  the 
precursor.  The  error  bars  denote  the  esti  mated  maximum  error  and  each 
plotted  point  is  an  average  of  from  three  to  five  separate  measurements. 


FIRINC  DELAY  <p>M> 

Fig.  1-36.  Normalized  impulse  delivered  to  0.38-cm-diameter  aluminum 
target  with  the  short  pulse  as  the  precursor.  The  error  bars  denote  the 
estimated  maximum  error  and  each  plotted  point  represents  an  average 
of  three  to  five  sc;»arate  measurements. 
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total  impulse  is  seen  to  be  just  beginning  to  recover  aftera  decay  of  1  msec.  Results  for  the 
case  in  which  the  short  pulse  was  the  precursor  are  shown  in  Fig.  1-36.  Here,  the  impulse 
delivered  to  the  aluminum  target  shows  a  steady  increase  and  takes  on  the  order  of  a  msec  to 
"recover." 

Previous  double  -pulse  impulse  measurements  on  small  lead  targets  presented  by  hall30 

indicated  about  a  50-percent  enhancement  of  the  impulse  delivered  by  the  second  pulse  after  about 

a  50-psec  delay.  Hall's  double-pulse  laser  system  was  such  that  each  pulse  had  an  energy  of 
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about  1  5  J,  pulse  lengths  of  about  25  psec,  and  power  densities  of  a  few  times  10  W/cm  incident 
on  the  target.  We  did  not  find  any  significant  enhancement  with  either  aluminum  or  carbon  targets 
although  it  should  be  noted  that  our  experiments  were  performed  with  power  densities  which  were 
an  order  of  magnitude  greater  and  less  t'  an  those  used  by  Hall. 

We  conclude  from  these  results  that  the  plasms  resulting  from  the  precursor  pulse  tends 
to  cause  the  energy  in  the  second  pulse  to  be  absorbed  in  a  large  volume  thereby  leading  to  a 
greatly  reduced  peak  pressure  at  the  target  and  also  c  “c-easing  the  total  momentum  transferred 
to  the  target.  This  view  is  supported  by  the  fact  that,  in  addition  to  the  reduced  peak  pressure, 
the  pressure  pulse  due  to  the  second  laser  pulse  was  both  delayed  and  broadened  in  the  double¬ 
pulse  mode.  The  delay  times  required  for  recovery  of  peak  pressure  and  impulse  to  those  ex- 
pec.ed  for  single-pulse  cases  were  fend  to  be  on  the  order  of  1  msec  for  the  experimental  con¬ 
ditions  reported  here. 

J.  E.  Lowder 
L.C.  Pettingill 


2.  Hydrodynamic  Effects  in  Double-Pulse  Measurements 

a.  Introduction 

To  supplement  the  pressure  and  impulse  measurements  of  the  resnonse  of  a  target  to  two 
10.6-pm  laser  pulses  (see  Sec.I-C-1)  we  have  undertaken  an  examination  of  the  hydrodynamic 
and  plasma  environment  through  which  the  second  pulse  must  travel.  We  have  further  examined 
the  velocity  of  the  luminous  front  of  the  second  pulse  and  correlated  its  motion  with  the  effects 
of  this  environment.  We  have  found  that  both  the  lowered  neutral  density  and  increased  electron 
density  in  the  region  of  the  target  lead  to  a  rapid  motion  of  the  luminous  front  away  from  the 
target  (LSD  wave  and  breakdown  wave).  This  rapid  motion  can  account  for  the  impulse  and 
pressure  decoupling  nnted  earlier. 

b.  Interferometry 

The  interferometer  used  ir_  this  experiment  is  the  same  Mach-Zehnder  interferometer  de- 
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scribed  in  the  Optics  Research  Report  with  only  the  5145-A  argon  laser  used  as  a  light  source. 
What  is  measured  then  is  the  contribution  to  the  index  shift  both  from  electrons  and  heavy  ions 
and  neutrals. 

Two  types  of  laser  pulses  were  examined.  The  first  pulse  in  -  double-puls**  experiment 
was  typically  about  40  psec  long,  at  an  intensity  of  1  to  5  X  108  W/cm3.  The  second  pulse, 
fired  some  time  delay  later,  was  about  1  psec  long  with  intensities  that  could  be  varied  up  to 
10^  W/cm2,  although  for  the  interferometric  investigation  intensities  of  5  x  107  to  2  x  108  W/cm2 
were  more  typical. 

Looking  at  the  long-pulse  data  first,  a  volume  heating  of  the  plasma  occurred  (see  Sec.  I-C-2d) 
generating  a  spherical  blast  wave.  Figure  l-37(a-b)  shows  interferograms  of  the  blast  wave  at 
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(a)  Two-dimensional  interferogram  (b)  Two-dimensional  interferogram 

for  1  =  1.3  MW/cm^,  At  =  27  jisec.  for  T  =  4.8  MW/cm^,  At  =  60  nsec. 


TIME  - ► 

(c)  Streak  interferogram  for  T  =  1.1  MW/cm2. 


Fig.  1-37.  Interferograme  of  laser-target  interaction  for  long-pulse  conditions. 
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two  different  intensities  and  times.  Figure  I-37(c)  is  a  streak  interferogram  for  the  conditions 
of  Fig.  I-37(a),  showing  the  index  growth  normal  to  the  target  as  a  function  of  time. 

Looking  at  the  short  pulse  alone  now,  Figs.  I-38(a-c.)  are  interferograms  at  different  times 
showing  the  evolution  of  an  initial  LSD  wave  [Fig.  I-38(a) )  into  a  cylindrical  blast  wave 
I  Fig.  1- 3 8 (b)  1  and  finally  a  sDherical  Hast  wave  [Fig.  1- 3 8(  c) ).  Note  that  at  late  times  the  shock 
detaches  from  the  hot  core  at  the  target. 

Finally,  Fig.  1-39  shows  an  interferogram  of  both  pulses,  with  a  time  delay  of  ~100  psec 
between  the  start  of  pulse  1  and  pulse  2,  where  the  interferogram  was  taken  110  psec  after  the 
start  of  pulse  1. 

An  Abel  inversion  program  has  been  written  to  reduce  the  interferometric  data  (see 
Sec.  l-C-4),  and  the  data  reduction  is  presently  under  way.  Figure  1-40  indicates  a  typical 
density  profile  for  the  long-pulse  data  taken  5  psec  after  the  beginning  of  the  long  pulse.  Notice 
the  density  at  the  target  is  negative  out  to  about  3.5  mm.  This  represents  the  electron  contribu¬ 
tion.  One  may  write  the  combined  contribution  of  the  free  electrons  and  neutral  density  to  the 
index  of  refraction  as: 

,,  n  l  Apt 

F  =  1.18  x  lo'22  +  2.77  x  10'4  -yD 

Here,  F  =  fringe  shift  in  wavelengths,  ng  =  electron  density,  l  =  path  length  of  electrons,  = 

ambient  density,  t  =  path  length  of  neutrals,  and  A  ~  interferometer  wavelength.  Using  this 

equation  and  assigning  the  entire  negative-density  contribution  to  the  electrons,  we  arrive  at  an 

electron -density  profile  and  nr  uiral-density  profile  as  shown  in  Figs.  J-41(a)  and  l-41(b),  re- 
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spectively.  There  appears  to  be  a  relatively  uniform  sphere  of  electrons,  with  ng  ~  4X10  cm  . 
This  is  consistent  with  the  two-wavelength  interferometric  measurements  previously  reported 
in  Kef.  3. 


c.  P  obe  Data 


To  examine  the  spatial  and  temporal  history  of  the  electron  distribution  in  front  of  ur.  alu¬ 
minum  target,  a  cross-beam,  10.6-pm  probe  laser  was  used  (as  shown  in  Fig.  1-42).  The  tem¬ 
porally  stable,  one-watt  CW  laser  beam  was  focused  to  0.1 -mm  diameter,  in  the  center  of  the 
plasma  column  which  is  typically  several  millimeters  in  diameter.  The  time-resolved  trans¬ 
mission  was  measured  1,  3,  5,  7.5,  10  and  15  mm  in  front  of  the  target  surface.  Figure  1-43 
shows  the  spatial  dependence  of  the  transmission  at  the  end  of  the  laser  pulse  and  20  psec  after 
the  end  of  the  pulse.  Figure  1-44  illustrates  the  temporal  recovery  of  the  transmission  3  mm 
in  front  of  the  target. 

If  one  makes  the  assumption  that  for  10 -pm  radiation  the  absorption  is  due  to  inverse  brems 
strahlung,  the  absorption  coefficient  can  be  used  to  determine  the  electron  density 


o(cm  *) 


-37  7  7- 

1.6  X  10  A  n  T 
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Here,  a  is  the  absorption  coefficient  in  cm  1 ,  A  is  the  probe  wavelength  in  pm,  n  is  the  electron 

_3  e 

density  in  cm  ,  Tg  is  the  electron  temperature  in  eV,  and  v  is  the  laser  frequency.  For 
10.6  pm  and  Tg  *  1  eV, 


n  *  2.5  x  10*  7  v/ 
e 
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(a)  LSD  initiation,  I  =  6  x  108  W/cm^,  (b)  Cylindrical  blast  wave  stage, 

At  =  1  nsec.  I  =  2.8  x  108  W/cm2>  At  =  40  nsec. 
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(c)  Spherical  blast  wave  stage,  T  =  2.2  x 
108  W/cra2,  At  =  iOO  nsec. 


Fig.  1-38.  Two-dimensional  intorferograms  of  laser-target  interaction 
for  short-pulse  conditions. 
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Fig.  1-44.  Plot  of  transmission  through  plasn.a  column  as  a  function 
of  time  for  2  =  3  mm.  The  filled  symbols  indicate  the  laser  is  still 
on.  The  open  symbols  are  appropriate  to  the  laser's  being  off. 

From  the  transmission  curves  of  Figs.  1-43  and  1-44,  we  can  derive  values  of  nl.  Assuming 
a  uniform  sphere  of  electrons  with  radius  of  about  i  cm  (i.e.,  where  the  absorption  is  -0  in 
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Fig.  1-43),  we  find  electron  densities  of  2  to  3  x  10"  /cm  at  the  target  for  times  from  40  tc 
60  psec  after  the  start  of  the  pulse.  The  numbers  are  consistent  with  the  values  found 
interfero  metrically. 

d.  Streak  Data 

A  series  of  experiments  was  run  to  determine  the  effects  of  the  neutral  and  electron  density 
left  by  the  first  pulse  on  the  dynamics  of  the  second  pulse.  For  these  experiments  the  first 
pulse  was  40  psec  long,  and  the  luminous  velocity  of  a  second  oulse  was  examined  as  a  function 
of  time  delay  between  the  two  pulses  and  intensity  of  the  second  pulse. 

Figure  1-45  shows  a  sequence  of  streak  photographs  at  time  delays  of  30  to  400  psec  between 
the  start  of  pulse  2  and  the  start  of  pulse  1.  For  the  early  time  history,  while  the  first  la»er 
pulse  is  still  on,  breakdown  apparently  initiates  at  the  edge  of  the  plasma  produced  by  the  first 
pulse,  and  a  normal  LSD  wave  is  launched  at  this  point.  In  addition,  enough  energy  is  trans¬ 
mitted  through  the  volume  absorbing  plasma  to  launch  a  weak  LSD  wave  from  the  target  surface. 
For  times  from  ~100  to  400  psec,  a  very  rapid  wave  is  launched  from  the  target  (v  ~  6  x 
106  cm/sec).  There  exist  two  possible  causes  for  this  high  velocity.  Firs*.,  it  might  be  an 
LSD  wave  launched  in  the  low-density  air  left  behind  the  first  pulse.  However,  since  the  velocity 
is  about  five  times  faster  than  might  be  expected  in  ambient  air,  this  would  indicate  a  residual 
density  of  ~0.008  amagat  (v  ~  l/p1/3).  It  is  hard  to  measure  the  density  to  this  accuracy  at  the 
targe*  because  of  the  electron  contribution  to  the  index  change. 

The  electrons  themselves  may  be  a  contributor  to  the  velocity  of  the  luminous  front  since, 
us  Smith31  has  shown,  the  threshold  for  breakdown  decreases  with  electron  density.  For  the 
high-intensity  second  pulse  (~10^  V'/cm),  a  breakdown  wave  can  be  launched  back  towards  the 
laser.  This  wave  will  continue  until  the  intensity  drops  to  the  point  where  it  is  below  the  break¬ 
down  threshc..’  I  for  that  distance  from  the  local  point.  At  this  point  a  normal  LSD  wave  will  be 
formed. 
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(a)  At  =  30  fisec.  (b)  At  =  1 50  nsec. 
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Fig.  1-45.  Sketches  of  streak  photographs  of  'luminous  pla?n?“u;/ron78 
resulting  from  double-pulsed  laser-target  interactions.  I}  ~  2  MW/cm  . 
T  ~  8  x  108  W/-.ir>.  Exposure  during  pulse  2  M  after  pulse  1. 


TIME  (fit #c) 


TIME  (fittc) 


(a)  I2  =  7.2  x  107  w/cm2 


(b)  I2  =  1.9  x  108  W/cm2 


(c)  I3  =  7.2  x  to8  W/cm2 


Fie  1-46  Sketches  of  streak  photographs  of  luminous  plasma  fronts 
resulting  from  double-pulsed  ’aser-target  Interaction-.  Exposure^unng 
second  pulse  that  occurs  85  psec  after  first  pulse.  Ij  ~  2  MW/cm  . 


I  ['Z/cm*  ) 

Some  evidence  as  to  which  mechanism  is  more  appropriate  can  be  obtained  fro.r  Fife.  1-4 <>, 

in  which  streak  data  are  shown  which  were  taken  at  a  fixed  time  (85  nsec)  after  the  sti  rt  of 

pulse  1  for  a  number  of  different  intensities.  Using  these  streak  data,  e  have  plotted  the  initial 

velocities  as  a  function  of  intensity  at  the  target  'Fig.  1-47).  Since  we  e:.pect  the  hydrodynamics 

to  be  the  same  for  all  three  intensities,  the  velocities  should  vary  as  I1/3  if  an  USD  wave  were 

launched  in  all  cases.  Notice  for  lower  intensities  the  l1^3  behavior  becomes  more  appropriate. 
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However,  for  intensities  above  3X10  W/cm  ,  the  dependence  on  intensity  is  much  stronger 
which  is  characteristic  of  a  breakdown  phenomena. 

The  LSD  velocities  in  the  lower  intensity  streak  photographs  indicate  a  residual  density  of 
0.2  to  0.5  amagat  85  psec  after  the  firing  of  pulse  1 .  Thus,  all  we  really  know  about  the  residual 
neutral  density  at  the  target  at  5  psec  is  that  the  density  is  0.2  amagat  or  less.  However,  the 
difference  between  assuming  vacuum  at  the  target  and  assuming  0.2  amagat  is  less  than  a  factor 
of  two  in  the  electron  density. 


e.  Conclusions 


We  have  examined  the  physics  involved  when  two  laser  pulses  interact  with  a  target.  In 
particular,  we  have  examined  the  target  environment  left  by  the  first  pulse.  What  is  left  at  a 
time  of  the  order  jf  10  to  100  psec  after  the  first  pulse  is  a  low-density  (<0.2  amagat)  electron- 
rich  (ne  08 1017/cm3)  environment.  For  second  pulses  with  intensities  of  ~108  W/cm2  an  LSD 
wave  is  launched  with  a  velocity  characteristic  of  the  low-density  environment.  For  intensities 
cf  10  W/cm  ,  a  breakdown  wave  is  launched  in  the  electron-rich  environment  at  the  lo'.v  break¬ 
down  threshold  characteristic  of  such  a  medium.  This  wave  travels  with  a  high  velocity  which, 
in  general,  will  be  a  function  of  geometry  of  the  experiment.  In  either  case,  there  is  a  strong 
decoupling  of  the  radiation  from  the  target  consistent  with  trie  pressure-impulse  results  reported 


elsewhere. 


H.  rileiman 
R.  W.  O'Neil 
L.C.  Pettingill 


3.  Empirical  Modeling  of  Vapor  Sphere  Expansion 
from  Laser-Exploded  Water  Droplets 

When  water  spheres  are  irradiated  with  a  pulsed  laser  beam,  the  energy  absorption  leads 
to  a  shock  formation  traveling  radially  out  from  the  drop  at  a  velocity  faster  than  the  sound  velocity 
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and  to  a  vapor  sphere  whit  h  follows  at  a  velocity  (initially  less  than  the  speed  of  sound)  which 
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monotonically  decreases  in  time  ’  If  a  precursor  laser  puts*  is  used  to  clear  a  fog,  a  sub¬ 
sequent  pulse  would  experience  beam  distortion  primarily  due  to  the  medium  index  fluctuations 
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induced  by  these  vapor  spheres.  Vve  have  used  the  data  assembled  by  Kafalas  et  al_.  ’  to 
parameterize  empirically  the  vapor  sphere  expansion. 

We  tested  as  possible  dependence  of  the  vapor  sphere  velocity  with  time,  functions  that 
yield  exponential  decays  and  decays  inverse  to  some  power  of  time  [i.e.,  v(t)  =  v^  e  '  or 
v(t)  =  v&(t  +  t0)  m  tj11,  where  v0,  At,  tQ,  and  m  are  parameters  to  be  determined!  We  then 
fitted  the  parameters  to  the  measured  values  cf  vapor  sphere  size  as  a  function  of  time  by  inte¬ 
grating  the  velocity  function  with  respect  to  time.  We  found  that  the  exponential  velocity  decay 
was  much  too  fast,  and  the  best  fit  for  the  inverse  time  decay  yields  an  optimum  fit  for  a  radius 
expansion  of  the  vapor  sphere: 

a(t)vapor  sphere  J  026  mm  ln(  +  rw 

where  t  is  in  psec  and  rw  is  the  initiri  water  vapor  size.  The  rate  of  expansion  as  measured 
from  the  experimental  data  display  -io  dependence  on  initial  droplet  size.  The  matching  has 
been  done  for  a  range  of  initial  particle  rad’i  and  for  times  out  to  about  40  psec  after  laser  pulse 
incidence. 

These  data  correspond  to  a  laser  pulse  of  10  MW/cm2  and  to  a  half-width  of  ~300  nsec. 
There  is  som  uTcation  that  the  time  dependence  scales  inversely  with  power,  i.e.,  as  the 
power  is  incre  sed,  ihe  rate  of  expansion  correspondingly  increases,  but  we  do  not  have  suffi¬ 
cient  data  to  make  a  quantitative  determination.  We  hope  soon  to  explain  the  observed  expansion 
rates  or.  a  theoretical  basis. 

H.  Granek 


4.  Abel  Inversion  Program 


In  support  of  the  interferometric  data  described  in  Sec.  I-C-2,  we  have  developed  an  Abel 
Inversion  Code  to  reduce  the  phase  shifts  to  density  contours.  For  a  ••ylindrically  symmetric 
density  variation  ir.  a  material  of  uniform  composition,  the  change  in  the  density  at  a  point  y 
away  from  the  axis  of  symmetry  is  given  by  Abel's  inversion  formula 
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where  fi  is  the  phase  shift. 

Here  it  is  assumed  the  change  in  the  density  is  related  to  the  change  in  the  index  cf  re¬ 
fraction  by 


where  and  PQ  are  the  ambient  index  and  density,  respectively.  The  data  from  interferograms 
(cf.  Sec.  l-C-2)  are  smoothed  by  a  cubic  spline  interpolation  and  *he  numerical  derivative  is 
subsequently  expanded  in  a  basis  set  of  linear  splines  defined  by 
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The  resulting  expansion  for  6'  is  given  by 


<5’(r)  =  £  w.(r)  6! 
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Substituting  into  the  inversion  formula  and  evaluating  the  remaining  derivatives  results  in  the 
following  expression  for  the  density  change  at  y 
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This  inversion  procedure  has  been  tested  with  Gaussian  phase  shift  variations  for  which  the 
inversion  integral  has  a  closed  form.  The  results  indicate  that  this  code  will  invert  phase  shifts 
with  an  error  limited  by  the  number  of  points  in  the  linear  spline  expansion  of  6’.  For  an  ex¬ 
pansion  containing  90  points  the  error  is  less  than  0.5  percent  -  well  below  that  needed  for  this 
analysis. 

Since  the  smoothed  data  a  e  similar  to  a  sum  of  Gaussian  curves,  it  is  assumed  the  resulting 
inversions  reflect  accurately  the  density  variations  for  such  data.  The  largest  sourc  of  nu¬ 
merical  error  in  the  whole  procedure  is  in  initially  smoothing  the  raw  data. 


D.  L.  Mooney 


D.  DEVICES 

1.  Dual  500-J  Electron-Beam  leaser 

The  dual  500-J  e-beam  laser  is  in  the  final  stages  of  checkout.  The  system  has  been  com¬ 
pletely  assembled  and  is  now  operational.  It  is  currently  utilizing  the  original  Lexan  and  the 
interim  Lucite  sustainer  cavities.  Aluminum  sustainer  cavities  have  been  delivered  and  will 
be  installed  early  in  1974. 

The  device  is  complete  except  for  the  automation  of  certain  auxiliary  sequencers.  A  new 
digital  fast  sequencer  to  provide  timing  between  the  two  laser  operations  and  to  facilitate  opera¬ 
tor  control  of  the  firing  sequence  is  ready  for  test.  Provisions  for  selectively  disabling  parts 
of  the  interlock  chain  for  ..est  purposes  will  be  provided  in  a  modification  now  being  designed. 

Each  laser  has  been  tested  independently  and  has  operated  successfully  with  good  beam 
quality.  Double-pulse  operation  has  been  achieved. 

The  Lincoln  Laboratory  two-pulse  laser  facility  has  been  constructed  to  study  the  propaga¬ 
tion  and  target  interaction  phenomena  of  multiple-pulse  lasers  for  an  arbitrary  irterpulse  sepa¬ 
ration.  For  common-path  problems  such  as  air  breakdown,  thermal  blooming  in  a  stationary 
air  column,  precursing,  fog  hole  boring,  and  subsequent  propagation,  two  250 -J  beams  are 
obtaineo  via  a  germanium  beam  splitter  (half  the  energy  is  lost  by  using  the  beam  splitter). 

For  applications  where  a  small  angular  offset  is  tolerable  (target  interaction)  or  desirable 
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(blooming  interaction  between  two  slewed  pulses),  two  500-J  pulses  are  available  as  indicated 
in  Fig.  1-48.  In  this  configuration,  beam  diagnostics  are  obtained  using  two  high-transmiitance 
gratings  angularly  displaced  to  provide  a  complete  independent  analysis  on  each  beam  In  the 
common-path  arrangement,  time-resolved  measurements  on  each  pulse  are  required.  The  op¬ 
tical  range  has  45  rr  cf  straight  path  or  90  m  from  a  single  reflection.  A  30-m-long,  50-cm- 
diameter  evacuable  range  is  avai.ible  for  controlled  atmospheric  propagation  or  vacuum  target 
interaction  experiments. 

A.  J.  Morency 

R.  W.  O'Neil 


2.  Other  Lasers 

a.  H^-HF  Electron-Beam  Laser  System 

The  objective  of  this  project  is  to  obtain  a  long-pulse  (1-  to  5-psec)  HF  laser  with  1  to 
10  joules  output  energy.  The  approach  will  be  to  excite  electrically  a  gas  mixture  containing 
H.,  and  HF  and,  therefore,  to  utilize  electron  pumping  of  with  a  subsequent  resonant  transfer 
to  HF  as  the  basic  inverting  process  of  the  HF  laser.  The  electrical  excitation  will  be  done  with 
a  cold-cathode  electron-beam-stabilized  discharge  system  of  commercial  manufacture. 

3 

The  basic  physical  principles  of  the  laser  were  described  in  the  last  QTS  and  in  Ref.  34. 

Here  we  will  describe  only  the  most  recent  progress  on  the  project  and,  in  particular,  work  on 
understanding  the  laser  using  both  a  theoretical  and  experimental  approach.  Also,  brief  mention 
will  be  made  of  the  preparatory  work  for  the  arrival  of  the  large  e-beam  system. 

(1)  Modeling  of  an  H2*HF  Discharge 

The  simplest  and  most  direct  way  to  excite  preferentially  the  H2  vibrational  levels  is  to  use 
an  electric  discharge  with  a  suitably  tuned  ratio  of  electric  field  to  number  density  (E/N).  For 
the  case  of  the  H.,-HF  mixture,  tunability  of  E/N  is  beat  obtained  with  an  e-beam  preionized 
sustainer  voltage.  Because  the  cross  section  for  vibrational  excitation  of  H^  by  electrons  is 
comparatively  small  (e.g.,  a  lactor  of  8  less  than  CO),  it  is  necessary  to  have  a  large  electron 
density  In  the  sustained  discharge.  Inspection  of  commercially  available  electron-beam  guns 
shows  that  a  cold-cathode  gun  gives  the  highest  integrated  current  density  'or  the  1  -  to  5-psec 
pulse  lengths  which  are  of  interest  to  Lincoln  Laboratory. 

Because  of  the  anticipated  high  rate  of  attacnment  in  HF  gas,  it  is  desirable  to  use  an  H^-HF 
mixture  which  contains  as  little  HF  as  possib’e  from  the  viewpoint  of  a  workable  small-signal 
gain.  Prior  work  on  kinetics  also  shows  that  such  a  low  HF  concentration  is  also  desirable  if 
one  is  to  have  a  slow  rate  of  vibrational-to-translational  energy  transfer  between  two  HF  mole¬ 
cules.  However,  use  of  such  a  low-Z  (atomic  number)  mixture  is  not  desirable  if  one  wishes  to 
achieve  a  high  ratio  of  electron  density  to  hydrogen  density,  a  i  Rio  which  is  necessary  for  rapid 
electron  pumping  of  the  H2  vibrational  levels. 

This  situation  may  be  rectified  in  our  case  by  addition  of  argon  (Z  =  19)  to  the  mixture.  Not 
only  will  argon  add  the  necessary  stopping  power  to  the  mixture,  but  it  also  will  leave  the  rates 
of  the  various  modes  of  vibrational  relaxation  essentially  unaltered.  This  results  from  the  fact 
that  an  inert  atomic  gas  can  relax  vibrational  levels  only  by  the  highly  nonresonant  process  of 
vibratlonal-to-translational  energy  transfer.  In  addition,  as  wi"’  be  discussed  in  the  following 
section,  argon  will  reduce  the  temperature  rise  in  the  H2  gas. 

In  order  to  estimate  the  electron  density  In  the  typical  gas  mixture  shown  in  Table  1-2,  it  is 
first  necessary  to  estimate  the  relative  rates  of  the  various  channels  of  electron  loss.  Diffusion 
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may  be  immediately  eliminated  because  of  the  high  total  gas  j.  ressure.  Because  of  their  low 
electronegativity  the  attachment  rates  of  argon  and  H2  are  vanishingly  small.  The  attachment 
of  HF  is  not  known;  however,  the  attachment  rates  in  HC1  and  HBr  have  been  measured  by 
Buchel'nikova.3^  If  we  conservatively  estimate  its  rate  to  be  the  same  as  HC1  and  we  use  the 
HF  gas  pressure  of  Table  1-2,  we  find  an  attachment  rate  of  1.4  x  10  sec  V 

The  recombination  rate  of  the  mixture  will  result  primarily  from  the  process 

e  +  A*  —  2A 

Its  characteristic  rate  will  be  (7  x  10“7)  (n  ).  If  we  assume  ng  »  0  x  to13  (to  be  verified  below) 
we  find  the  characteristic  rate  to  be  3  x  to7  sec-1.  It  is  apparent,  therefore,  that  recombination 
is  the  primary  means  of  electron  loss  for  the  discharge  under  consideration. 

A  consistency  check  of  the  above  discussion  can  now  be  made  by  verifying  that  an  n  ~ 

13  e 

5x10  is  a  reasonable  electron  density  to  attain  with  the  electron-beam  and  discharge  param¬ 
eters  envisaged  here.  Since  the  recombination  time  is  fa^t  compared  to  the  inverse  of  a  typical 
pulse  duration  time,  we  may  assume  a  steady-state  electron  density  exists;  then,  ng  »  vS/o 
where  S  is  the  electron  production  rate  and  or  is  the  recombination  coefficient.  Using  the  tabu¬ 
lated  stopping  power  of  argon,  we  car.  calculate  S  «  lO^1  electrons/second.  For  argon  molecuiar 
recombination,  a  =  7  x  10  7,  and  thus  ng  *  5  x  to13  cm3  is  reasonable. 

As  shown  in  Ref.  34,  the  HF  gain  is  sensitive  to  both  the  vibrational  and  translational  rota¬ 
tional  temperature  of  the  gas  mixture.  Further,  the  rate  and  fractional  power  transfer  to  each 
of  these  modes  depends  nonlinearly  on  the  E/N  of  the  discharge.  As  a  result,  it  is  not  possible 
to  calculate  the  magnitude  or  temporal  behavior  of  gain  without  a  complete  solution  of  the  time- 
varying  problem.  Here  we  describe  the  results  cf  a  computer  modeling  of  the  discharge  which, 
while  not  yet  complete,  is  sufficiently  developed  to  allow  Dinpointing  of  the  specific  problems  in 
the  laser  design  as  well  as  their  possible  solutions. 

The  basic  element  of  our  computation  is  calculation  of  the  electron  energy  distribution  in 

the  discharge.  The  overlap  of  this  d’s.ribution  with  the  cross  sections  of  each  of  the  processes 

of  interest  will  determine  the  rate  of  that  process.  To  determine  the  distiibution  we  have  solved 

the  Boltzmann  equation  much  in  the  manner  of  Frost  and  Phelps3^  (and  later  Nighan37).  The 
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vibrational  cross  sections  used  in  this  calculation  are  those  of  Schultz;  all  other  cross  sections 
are  from  Frost  and  Phelps.  Except  where  noted,  the  physical  parameters  used  i-i  the  calculation 
are  those  given  earlier  in  Table  1-2. 

Figure  1-49  shows  the  initial  rate  of  electron  pumping  of  the  H-,  vibrational  levels,  i.e.,  the 

pumping  rate  assuming  no  excited  H2  molecules,  for  various  values  of  E/N.  From  this  figure 

it  is  apparent  that  if  (me  wishes  to  pump  on  the  scale  of  10  psec  or  less,  then  an  E/N  of  less  than 
-16  2 

10  V/cm  is  unacceptable. 

However,  in  order  to  obtain  an  accurate  representation  of  the  pumping  rate,  it  is  essential 
to  take  into  account  the  change  in  the  electron  distribution  with  the  degree  of  excitation  in  the 
gas.  This  arises  because  of  the  fact  that  as  the  population  in  the  first  vibrational  lev»l  rises, 
superelastic  collisions  can  then  give  energy  back  to  the  electrons. 

Thus,  to  determine  how  the  vibrational  population  changes  in  time,  we  must  solve  a  nonlirsar 
rate  equation  such  as 

dH,(v  =  1)  r  r 

- - =  H2(v  =  0)  nfi  \  f(u,  Nq,  Nt)  <r(u)  udu  -  H2(v  =  1)  n&  \  f(u.  Np.  Nt)  o'(u)  udu  (1-6) 
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Fig.  1-49.  Rate  of  electron  pumping 
of  the  H2  vibrational  levels  assuming 
an  H?  vibrational  temperature  of 
300"K  as  a  function  of  electric  field 
to  neutral  particle  density. 


Fig.  1-50.  Calculated  time -dependent 
gain  of  the  Pj(8)  transition  assuming 
0.1  torr  of  HF  and  200  torre  of  H2  as 
a  function  of  electric  field  to  neutral 
particle  density. 


where  Nn  .  is  population  in  the  zeroth  and  first  vibrational  levels,  respectively;  u  is  the  electron 
0,1  l/2 

energy  in  Ev;  (u,  Ng,  Nj),  (u)  '  is  the  electron  energy  distribution  (a  function  of  u  and  vibra¬ 
tional  level  population);  and  a(u)  and  <r'(u)  are  the  cross  sect<ons  for  electron  excitation  an;  de¬ 
excitation  of  the  vibrational  levels,  respectively.  If  the  degree  of  vibrational  excitation  is 
small,  we  may  neglect  all  excited  levels  higher  than  the  first,  and  (1-6)  is  the  only  rate  equation 
to  be  solved. 

Notice  that  since  the  vibrational  equilibration  time  for  H2  and  HF  is  rapid  (300  nsec)  com¬ 
pared  to  tiie  anticipated  pumping  pulse  length,  we  can  determine  the  population  of  HF  in  the  first 
vibrational  level  by  assuming  intermolecular  vibrational  equilibr.um  exists.  The  results  of  a 
calculation  of  HF  gain  on  the  P(8)  v  =  t  ■»  v  =  0  transition  as  a  function  of  time  for  three  values 
of  E/N  are  shown  in  Fig.  1-50.  It  is  apparent,  first  of  all.  that  the  rate  of  excitation  is  an  ex¬ 
tremely  rensitive  function  of  E/N.  Therefore,  it  is  vital  to  have  the  freedom  to  tune  the  dis¬ 
charge  E/N  as  would  exist  in  an  e-beam  sustained  discharge.  Second,  the  limitation  on  low 
values  of  E/N  is  not  simply  pumping  time.  Rather,  the  fraction  of  discharge  energy  which  goes 
into  gas  heating  (i.e.,  translational  and  rotational  energy)  increases  as  one  goes  toward  lower 
values  of  E/N.  As  a  result,  the  tendency  is  for  a  higher  gas  temperature  (and  hence  lower  gain) 
for  a  given  vibrational  temperature  as  one  decreases  E/N.  This  is  illustrated  in  Fig.  1-50  by 
the  smaller  gain  obtained  in  the  curve  for  E/N  =  1.5  X  10  ^  V/cm2  than  that  for  2  x  10  1  ^  V/c n2. 
A  third  observation  is  that  if  one  goes  toward  higher  E/N,  discharge  energy  is  then  wasted  on 
the  next  higher  inelastic  process  (HF  and  H2  dissociation). 

The  problem  of  gas  heating  is,  in  reality,  a  more  serious  problem  than  shown  in  Fig.  1-50. 
Our  model  shows  that  in  a  pure  H2-HF  mixture  all  P-branch  transitions  with  J  <  8  do  not  exhibit 
gain  primarily  because  of  the  effects  of  gas  heating.  Hence,  control  of  gps  heating  is  considered 
a  major  experimental  difficulty  in  design  of  this  laser.  If,  however,  one  allows  for  the  presence 
of  argon  in  the  gas  mixture,  the  problem  of  gas  heating  is  3een  to  be  much  less  severe.  Since 
the  first  inelastic  cross  section  in  argon  occurs  at  11  eV  and  since  tne  average  electron  energy 
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EXPEP  ENT  WITH  ELECTRONS  PHOTOIONIZED 


Fig.  1-52.  Apparatus  used  for  scaled-down  modeling  of  H^-HF  laser  system. 
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for  a  typical  H..-HF  plasma  considered  above  is  1  eV,  at,;-.  can  effect  such  a  discharge  only 

via  its  elastic  cross  section  (or  cross  section  for  momentum  ransfer).  However,  in  the  region 

of  1  eV  or  less,  the  Ramsauer  minimum  in  the  argon  elastic  cross  section  makes  it  much  smaller 

than  the  corresponding  cross  section  in  As  a  result,  for  a  gas  mixture  of  equal  parts  of  H 

and  A  for  an  E/N  ratio  giving  an  average  energy  of  1  eV  or  less,  the  presence  of  A  has  no  effect 

on  the  electron  energy  distribution.  Further,  there  ic  no  heating  of  argon  by  momentum  transfer 

collisions  with  electrons  because  M  /M  «  M  /Mu  where  M  is  the  mass  of  the  electron,  etc. 

e  a  e  “2  e 

Thus,  the  chief  effect  of  adding  argon  is  to  increase  the  heat  capacity  of  the  gas  mixture.  The 
dramatic  effects  of  adding  argon  to  the  gas  mixture  are  shown  in  Fig.  1-51.  We  see  that  by  de¬ 
creasing  the  rate  of  temperature  rise,  the  addition  of  arf  on  causes  the  gain  to  rise  faster  and 
to  attain  a  greater  peak  magnitude  t  ian  without  it. 

In  conclusion,  we  again  stress  tnat  this  study  is  by  no  means  complete  yet.  To  be  included 
in  our  model  art  a  more  detailed  analysis  of  the  effects  of  an  argon  additive,  recently  published 
revised  cross-section  data  for  H2  vibrational  and  rotational  levels,  and  investigation  of  the  ef¬ 
fects  of  considering  more  than  two  vibrational  levels  for  both  the  H2  and  HF  vibrational  manifold. 
We  summarize  the  results  of  this  section  as  follows: 

(a)  Witi...i  the  limitations  of  our  present  model,  small-signal  gains  of  practical 
magnitude  (~0. 002/cm)  can  be  obtained  in  a  cold-cathode  electron-beam 
sustained  discharge  containing  H2<  HF,  and  argon. 

(b)  The  most  serious  difficulty  which  has  become  apparent  in  the  study  thus 
far  is  that  of  gas  hcat;ng.  However,  its  effects  can  be  reduced  by  judicious 
addition  of  argon. 

R.  Osgood 
D.  L.  Mooney 

(2)  Experimental 

Although  every  attempt  has  been  made  to  include  all  significant  detail  in  the  above  numerical 
calculations,  uncertainties  as  to  the  particular  values  of  many  cross  sections  make  it  advisable 
to  check  the  predictions  of  the  model  with  a  small-scale  experimental  apparatus.  Such  an  ap¬ 
paratus  has  heen  constructed  and  is  shown  schematically  in  Fig.  1-52.  Note,  first,  that  the  system 
is  constructed  of  components  resistent  to  chemical  reaction  by  HF.  Gases  are  flowed  through 
the  system  so  that  ar.y  discharge  byproducts  will  not  remain  in  the  discharge  region  after  pne 
pulse  cycle.  A  cold  trap  is  used  to  remove  HF  from  the  H2”HF  mixture  prior  to  its  entering 
the  vacuum  pump. 

Electrical  excitation  of  the  gas  mixture  is  accomplished  by  applying  a  DC  sustainer  voltage 
to  the  two  stainless  steel  electrode  ^tates  within  the  sample  cell  and  then  generating  a  pulse  of 
electrons  within  the  gas,  thus  a]  owing  current  to  flow.  Since  one  of  the  purposes  of  this  appara¬ 
tus  is  .o  simulate  a  larger  cold-  cathode  system,  the  ideal  way  tc  generate  the  electrons  within 
the  g^s  is  to  utilize  a  small  cold  -cathode  e-beam.  Currently,  such  a  device  is  being  tested. 

If  it  works  satisfactorily,  it  will  ie  incorporated  in  the  above  apparatus.  Meanwhile,  preliminary 
experimentation  has  been  done  using  photoionization  with  a  tri-n-propylamine  seeded  mixture 
While  such  a  seeded  mixture  is  not  ideal  from  the  standpoint  of  chemical  stability,  it  can  provide 
useful  information  on  rates  of  attachment  ar  recombination  in  H.,-HF  mixtures.  As  mentioned 
in  the  previoua  section,  such  information  is  necessary  if  one  wishes  to  calculate  the  electron 
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density  in  an  e-beam  preionized  plasma.  The  results  of  the  experiments  with  the  photoionized 
mixture  show  that  for  the  typical  gas  mixtures  given  in  Table  1-2  dissociative  recombination 
with  argon  molecular  ions  is  the  process  limiting  the  electron  density.  Electron  attachment  to 
HF  does  not  appear  to  be  a  competitive  process. 

R.  M.  Osgood 

(3)  Preparation  for  Large  Cold-Cathode  System 

To  prepare  for  the  arrival  of  the  large  e-beam  system  from  the  subcontractor,  a  number 
of  pieces  of  auxiliary  apparatus  have  been  designed,  fabricated,  and  tested.  Among  these  are 
two  transverse-discharge  HF  lasers.  These  will  be  used  as  probes  for  measurement  of  gain 
in  the  larger  high-pressure  device.  Also  fabricated  were  various  components  of  the  gas-handling 
and  vacuum- exhaust  apparatus  for  that  device.  Much  of  the  experience  gained  in  building  the 
gas  flow  sv.iem  for  the  smaller  ^‘HF  system  (described  above)  has  been  of  use  in  this  phase. 

R.  M.  Osgood 

A. J.  Morency 


b.  E-Beam-Initiated  HF  Laser 

The  Laser  Technology  program  requires  hign-energy  pulsed  HF  lasers  to  perform  beam- 
propagation  and  surface-effects  experiments  in  the  3-jim  region  of  the  spectrum.  In  order  to 
determine  the  effectiveness  of  a  repetitively  pulsed  laser  within  the  limitations  of  our  budget, 
these  experiments  will  be  performed  with  two  pulses,  from  two  distinct  lasers,  of  adjustable 
time  separation.  One  group  of  experiments  requires  two  pulses  of  approximately  200  joules  and 
200 -nsec  length.  We  plan  to  achieve  this  goal  by  using  high -intensity,  high-energy  electron 
beams  to  initiate  HF  chemical  lasers. 

We  are  currently  building  the  first  of  these  lasers.  The  plan  is  to  introduce  a  500 -keV, 
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60-nsec  electron  beam  of  20  cm  cross  section  and  250  kA  peak  current  at  a  small  angle  to  the 
laser  axis.  The  beam,  which  will  be  bent  and  confined  by  an  axial  magnetic  field,  will  initiate 
the  reaction  by  dissociating  a  fluorine -bearing  compound.  Initially,  the  laser  mixture  will  be 
SF^  and  C,H^  from  which,  based  on  work  reported  at  Maxwell  Laboratories  and  Sandia  Labo¬ 
ratories,  0  we  can  expect  at  least  60  joules.  We  have  ordered  the  electron-beam  accelerator, 
which  will  be  delivereu  in  March  1974.  The  laser  and  associated  components  are  currently 
being  fabricated  and  will  be  ready  by  that  time.  Unstable  resonator  geometry,  which  has  been 
shown  to  be  useful  for  high-gain,  large-volume  lasers,  has  been  incorporated  into  the  design. 

It  is  expected  that  the  laser  will  be  operational  within  two  months  subsequent  to  delivery  of  the 
electron-beam  accelerator.  To  ensure  maximum  field  homogeneity,  a  computer  program  has 
been  written  to  determine  the  optimum  nu'.iber  of  windings  per  individual  coil. 

1 

S.  Marcus  R.C.  Archibald 

A.J.  Morency  D.  L.  Mooney 

W.  R.  Fanning 

c.  HF  Electric-Discharge  Laser 

Work  was  continued  this  period  on  the  eventual  development  of  an  electric -discharge  HF 
laser  of  10  to  20  joules  near  diffraction -limited  output.  Due  to  excessive  arcing  in  the  discharge, 
the  maximum  output  obtained  until  this  period  was  0.13  joule  from  a  0  2 -liter  dev-ice  with 
100  joules  input  -  less  than  10  percent  of  that  expected.  Efforts  this  period  were,  therefore, 
concentrated  on  solving  the  arcing  problem. 

\ 


\ 
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It  was  decided  to  use  smooth  meta'J’c  e'.-trodes,  shaped  for  maximum  field  uniformity, 
and  to  improve  the  discharge  hy  (1)  minimizing  the  combined  inductance  of  the  laser  and  current 
source,  and  (2)  adjusting  the  laser  mixture  so  that  its  impedance  matches  that  of  the  current 
source,  thus  minimizing  the  time  required  for  energy  transfer  to  the  gas  and  reducing  the  proba¬ 
bility  of  arcing. 

The  most  dramatic  improvement  in  laser  periorma.ice  resulted  from  replacing  as  a 
hydrogen  atom  source  with  C^H^.  V-'m  a  mixture  of  typically  65  torrs  SF^  and  10  torrs 
and  an  electrical  input  of  100  joules,  an  output  of  2.25  joules  was  obtained.  This  efficiency  is 
comparable  to  the  best  obtained  from  nonchain  HF  lasers. 

During  the  next  period,  we  plan  to  test  a  commercially  obtained,  low-inductance  Marx  bank 
with  >500  joules  storage  in  an  attempt  to  increase  further  the  output  of  the  laser.  We  also  plan 
to  use  an  unstable  resonator  to  optimize  its  beam  quality. 

S.  Marcus 


d.  H2-F2  Laser 

It  is  possible  that  the  mixture  of  SF^  and  C^H^  will  not  yield  sufficient  energy  to  provide 
200-joule  laser  pulses  which  are  needed  for  3-pm  surface  effects  and  propagation  experiments. 

In  that  event,  we  will  be  compelled  to  use  the  H2*F2  chain  reaction,  from  which  more  energy 
is  available.  With  this  reaction,  laser  pulses  cam  be  stretched  at  a  cost  of  output  energy  to 
cover  the  1-  to  5-psec  pulse-length  regime  oy  adjustment  of  gas  pressure.  Because  of  this, 
such  a  laser  can  also  serve  as  a  backup  for  the  H2*HF  laser  described  in  Sec.  I-D-2a  above. 

Use  of  such  a  corrosive  and  explosive  medium,  however,  will  require  stringent  construction 
techniques  and  safety  procedures.  In  order  to  develop  an  expertise  for  handling  H2"F2  mixtures, 
we  plan  to  build  a  small  flishlamp-initiated  1L,-F2  laser  of  roughly  10  joules  output  capability. 

To  date,  we  have  received  almost  all  of  the  required  components  for  such  a  system,  and  we 
shall  begin  assembly  in  the  near  future. 

S.  Marcus 
H.  Granek 


e.  HF  Laaer  Kinetics 

In  support  of  the  high-energy  HF  laser  device  program,  a  computer  study  of  the  molecular 
kineticc  of  the  laser  has  been  initiated.  The  goal  is  to  be  able  to  predict  the  laser  output  as  a 
function  of  the  relevant,  parameters.  The  following  assumptions  constitute  the  model  we  have 
chosen: 

(1)  F  atcms  are  produced  by 

SF^  +  e  —  SFg  +  F  +  e 

(2)  Vibrationelly  excited  HF  is  produced  by 

F  +  H2  -  HF(v)  +  H  . 

(3)  Only  HF  levels  with  v  <  3  are  produced  and  relaxation  processes  involving 
v  >  3  are  negligible. 
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(4)  The  dominant  relaxation  processes  are: 

HF(v)  +  HF(v')  ^  HF(v  H-  1)  +  HF(v'  -  1) 
HF(v)  +  HF  -*  HF(  v  -  1)  +  HF 
HF(v)  +  F  -  HF(v  -  1)  +  F 
HF(v)  +  H2(0)  -  HF(v  -  1)  +  H,(l) 

(H2(0)»H2) 


(5)  Rotational  relaxation  is  sufficiently  fast  so  that  the  rotational  energy 
distribution  is  thermal  with  the  rotational  temperature  equal  to  the 
translational  temperature. 

Under  these  assumptions,  the  controlling  differential  equations  are: 

(D  ^  =  -f  J<t)  (SF6)  -ktH2F 
d(SF,) 

(2)  ___<L  =  -  f  J(t)  (SF6) 
dH? 

<3>  ~dT  -  -W 

(4)  ~3  -  k3H2F  -  X3n3n  -  y^F  +  k.2n2  -  k^n,  +  k^n,^ 

-  kbl2nQn3  -  kH^n3H2  -  S32 


dn2  2  2 
(5)  =  k2H2F  +  X3n3r,  -  *2n2n  +  Y^n^  ~  Y2n2F  +  kfinj  “  2kf2n2 

+  2kb2n3n3  ~  kfi 2n* H  kbl2n0n3  +  kH23n3H2 

“  kH2zn2H2  +  S32  _  S21 


dn 


1  2  2 
(6)  — —  =  kjH-jF  +  X2n2n  -  Xj^n  +  y2n2F  -  YjHjF  +  kf2n2  -  2k^nj 


+  2kbln0n2  k-b2nln3  kfl  2ni  a2  +  kbl2n0r‘3  +  kH22n2H2 


k|3  n,  H,  +  S„ ,  S, 


21 


12  21  10 


c"n0  2 
(7)  dT  =  k0H2F  +  Vln  4  r^jF  +  kfln,  -  k^m,  +  k^n^ 


"  ltbl2n0n3  +kH21nlH2  *  S10 


where 

n  =  n^  +  n3  +  n2  +  n3 
k^  =  kj  +  kz  +  kj  + 
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and 


wheri'  J  +  J 

v-1  v 


(8),  (9),  (10) 


Also, 


Definitions: 


t 

F 

<T 

J(t) 


e 


n 


r, 

*fv 


V,V‘1  4?kTE 


\  .  *  .  B!2J  +  i) 

v,v-1  v,v-l  v 


,ir  ,  t 
V.v-l  v.v-l  V, V - 1 


/- BJV-1(JV-1  +  1>\ 

-nv-i  e*P\ - kT - j 


nv  exp 


PM 


1  1  for  v  =  1,  2,  and  3  are. 


d* 


v,v-l  . 

-n -  =  c<t  , 

dt  v,v-l 


t\  .  B(2  J  +  1) 
V,  v-1  V 


nv  expl 


+  id  ln<‘  "K> 


,  T 

.kT 

V, 

/-I  V 

.v-1 

r  (J 

V  V 

+  lh 

-  n  . 

v-t 

kT 

1 

an 

E 

+  f/ 

V 

v.v-l 

dt 

T  . 

v,  v-1 

pj;  (2Jv  +  1)  exp 


/-BJ  (J 

I  v  v 

\  kT 


*  =  *32  +  *21  +  *10 


out  2  -  K 


time 

fluorine  atom  concentration 

cross  section  for  dissociation  of  F  from  SF^  by  an  electron 

electron  current  density 
electron  charge 

total  rate  constant  for  F  +  —  HF  +  H 

concentration 
SF^  concentration 

concentration  of  HF  with  vibrational  quantum  number  v 
rate  constant  for  F  +  H^  —  HF(v)  -i  H 
total  HF  concentration 

rate  constant  for  HF'v)  +  HF  -»  HF(v  -  1)  +  HF 

rate  constant  for  HF(v)  +  F  -»  HF(v  -  1)  t  F 

forward  rate  constant  for  2  HF(v)  Z  HF(v  -  1)  +  HF(v  +  1) 
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*Sv 

kfl2 


*^12 

X 

V,  v-1 


B 

Jv 

k 

T 

E 


v.v-t 


Av 


v,v-l 


backyard  rate  constant  for  2  HF(v)  Z  HF(v  -  1)  +  HF(v  +  1) 
rate  constant  for  H2(0)  +  HF(v)  -  H2(l)  +  HF(v  -  1) 

forward  rate  constant  for  HF(i)  +  HF(2)  -  HF(0)  +  HF(3) 

backward  rate  constant  for  above 

wavelength  of  v  -*  v  -  1  transition 

laser  intensit;  inside  cavity  in  v  -»  v  -  1  transition; 

*  =  total  intensity  in  cavity;  *QUi  =  output  intensity. 

rotational  energy  constant 

rotational  quantum  number  of  upper  state  of  v  —  v  -  1  transition 

Boltzmann  constant 
translational  temperature 
energy  of  v  -»  v  -  1  transition 

Lorentz  width  of  v  —  v  -  1  transition 


Tv  v_j  spontaneous  emission  time  constant  for  v  —  v  -  1  transition 

c  speed  of  light 

f  length  of  gain  medium 

d  cavity  length 

K  output  coupling  fraction 

<*>e  vibration  energy  constant  for  HF 

<*>exe  anharmonic  correction  to  HF  vibrational  energy 

f  fraction  of  spontaneous  emission  in  solid  angle  contributing 

to  oscillation  buildup 

A  program  is  presently  bei'ig  written  to  solve  this  set  of  first-order,  second-degree  differ¬ 
ential  equations.  The  Runge-Kutta  integration  method  is  being  used. 

S.  Marcus 
D.  L.  Mooney 


f.  C02  Isotope  Laser  Measurements 

We  have  previously  reported  a  brief  summary  of  our  experiments  on  12C180,,  13Cl60,, 

1318  LL 

end  C  02  isotope  lasei  Detailed  description  of  these  experiments  and  the  resu'ts  will  be 

published  in  the  February  1974  issue  of  the  Journal  of  Molecular  Spectroscopy. 

We  have  obtained  laser  action  in  *2C^80*80  and  ^3C*80*80  isotopes  and  obtained  some 
preliminary  data.  Because  of  the  vast  number  of  lines  obtainable  In  the  mixed  isotopes,  identi¬ 
fication  is  a  more  iedibus  process.  To  alleviate  this  problem  and  «>  obtain  more  accurate 
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measurements,  the  final  results  will  be  obtained  with  the  lasers  locked  to  the  standing-wave 
saturation  resonances  observed  in  low-pressure  CO£  absorption  cells,  utilizing  the  upper  to 
ground  state  spontaneous  emission  band  at  4.3  pm.  Appropriate  external  cell  stabilization  set¬ 
ups  have  been  designed  and  built. 

The  design  of  the  grating-controlled  lasers  has  been  improved.  The  previously  troublesome 
output  mode  distortions  have  been  eliminated,  and  also  a  large  numoer  of  additional  lasing  tran¬ 
sitions  were  observed. 

C.  Freed 
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li.  OPTICAL  MEASUREMENTS  AND  INSTRUMENTATION 


A.  TRIANGULAR  SCANNING  MTF  INTERFEROMETER 

1.  MTF  Interferometric  Measurements 

Twc  engineered  models  of  the  triangular  scanning  MTF  interferometer1  have  now  been  built 
to  our  design  by  Infrared  Industries  for  M.l.T.  Lincoln  Laboratory.  It  is  suitable  for  field  use, 
with  (.-ovision  for  simple  visual  alignment  by  means  of  a  reticle  in  an  eyepiece.  A  camera  can 
be  attached  to  photograph  interferograms.  A  dove  prism  to  rotate  the  wavefront  has  been  fitted 
in  front  of  the  entrance  aperture.  A  flip  mirror  switches  from  visual  examination  to  photoelec¬ 
tric  detection.  The  shear  plate  in  normal  operation  is  continuously  rotated  at  3600  rpm  by  a 
motor,  but  can  also  be  manually  tcrred  or  locked  in  position  at  any  particular  shear  value  (motor 
switched  off).  An  optical  encoder  Keeps  track  of  the  angular  position  of  the  shear  plate,  provid¬ 
ing  indexing  pulses  once  per  revolution,  together  with  a  continuous  200- kHz  digital  output.  These 
two  encoder  data  allow  indexing  of  the  zero  shear  position  in  the  MTF  data  as  well  as  very  accu¬ 
rate  shear  (spatial  frequency)  axis  calibration. 

The  MTF  signal  data  and  the  two  encoder  data  outputs  are  recorded  on  magnetic  tape  (having 
a  high-frequency  bandwidth),  and  the  tapes  subsequently  digitized  and  fed  to  a  computer  for  data 
reduction  and  statistical  analysis  of  the  data.  The  instrument  measures  and  records  14,400  MTF 
curves  each  minute,  and  a  complete  MTF  curve  is  made  i.i  1.3  msec.  A  micrometer  adjustment 
allows  accurate  rotation  of  the  optical  part  of  the  triangular  interferometer  as  a  unit  such  that 
the  carrier  signal  can  be  easily  adjusted  while  operating  the  interferometer  (from  zero  to  about 
100  kHz).  The  instrument  has  initially  been  designed  for  use  with  a  HeNe  CW  laser  source 
(632.8  nm).  However,  within  certain  spectral  and  spatial  coherence  limitations,  it  has  been 
successfully  tested  with  a  small  white  incoherent  source,  using  suitable  bandwidth  filters.  We 
therefore  expect  to  be  able  to  look  at  star  sources  and  measure  the  MTF  at  various  slant  paths 
through  the  whole  of  the  atmosphere.  The  instrument  can  be  converted  lo  infrared  wavelengths, 
and  this  is  being  done  for  a  1  C-micron  version. 

A  photograph  of  the  instrument  is  shown  in  ^ig.  11-1.  Modular  amplifiers  and  filters  are  used 
to  provide  simple  electronic  processing  (mainly  noise  reduction  and  gain  control)  of  the  data  sig¬ 
nals  for  recording  on  tape.  The  interferometer  accepts  up  to  a  10-mm-diameter  collimated  light 
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beam.  As  an  illustration  of  the  MTF  output  from  the  interferometer,  a-1  89-mm-diameter 
Questar  telescope  was  used  as  a  receiver,  and  a  laser  beam  propagated  over  a  14-meter  path 
in  the  air-conditioned  laboratory.  Figure  11- 2(a)  shows  the  type  of  output  obtained  recorded  on 
tape  (6  MTF  curves  are  shown),  and  for  detailed  examination,  can  be  expanded  as  shown  Ly  the 
series  in  Fig.  ll-2(b).  Typically,  fluctuations  in  the  temperature  in  the  laboratory  air  produce 
the  variations  shown  between  the  curves  (these  three  pairs  of  curves  in  Fig.  11- 2(b )  were  selected 
at  about  one-minute  intervals).  The  curves  change  slowly  under  these  conditions.  If  now  labora¬ 
tory  air  is  blown  across  in  front  of  the  aperture  of  t  ie  telescope,  the  MTF  curve  does  not  signif¬ 
icantly  change  its  characteristics,  as  shown  in  Fig.  11- 3(a).  However,  if  temperature  fluctuations 
are  now  introduced  into  this  layer  of  air  blown  in  front  of  the  telescope  (by  including  a  heated 
coil  in  the  blower),  the  MTF  curve  collapses,  as  shown  in  the  three  pairs  of  curves  of  Fig.  11-  3(b). 
The  three  pairs  of  curves  were  selected  at  about  one-minute  intervals. 

In  summary,  large  masses  of  slow-movi.ig  air  from  the  air  conditioner  and  people  moving 
about  the  laboratory  are  causing  slow  variations,  but  not  drastic  collapse,  of  the  MTF  curve  in 
Fig.  11-  2(b).  A  fast  flow  of  ambient  air  (from  a  hair  dryer  located  three  feet  away),  about  two 
feet  in  front  of  the  telescope,  produces  no  change  in  the  nature  or  level  of  the  MTF  curve  unless 
the  heat  is  turned  on  to  introduce  large  rapid  temperature  fluctuations,  as  shown  in  Fig.  11- 3(b). 

In  this  latter  case,  the  MTF  collapses  and  fluctuates  (at  high  shear)  rapidly.  The  two  main  ex¬ 
amples  illustrated  are  intended  to  demonstrate  the  type  of  effect  on  the  MTF  produced  by  long- 
path  propagation  through  the  atmosphere  and  propagation  through  a  flight  boundary  layer,  respec¬ 
tively.^  Clearly,  a  more  detailed  simulation  and  analysis  is  called  for  in  order  to  obtain  realistic 
and  accurate  results.  We  are  using  statistical  methods  of  analysis  applied  to  the  large  number 
of  MTF  curves  that  can  be  measured  in  this  way,  to  evaluate  "mean,"  "peak,"  "minimum"  MTF 
curves,  and  variance  of  the  MTF  amplitude  at  each  shear  value  in  order  to  characterize  both 
the  boundary  layer  in  the  flight  tests  and  the  long-path  propagation  problem. 

2.  Angle-of- Arrival  and  Nonimaging  Tracker  Application  of  the  Interferometer 

Optical  vrackers  are  usually  based  upon  controlling  and  measuring  the  position  of  an  image 
on  a  suitable  sensor  by  some  optical  means  (such  as  mosaic  detector  arrays  or  rotating  reticles). 
However,  an  alternative  interferometric  nonimaging  technique  can  be  used  and  will  be  outlined 

here.  The  angle  of  arrival  of  an  optical  wave  can  be  measured  or  tracked  by  means  of  the  differ- 

3  4 

ential  path  change  in  a  two- beam  shearing  inte'-fe»,ometer.  (For  example,  see  Baker  and  Dyson.  ) 

In  practice,  the  corner  cube5  and  the  triangular1  shearing  interferometers  can  be  used  for 
angle-of-arrival  measurements  and  for  tracking  objects.  The  pointing  sensitivity  of  the  corner- 
cube  interferometer  has  been  analyzed  by  the  writer5  and  particular  application  of  this  technique 
to  tracking  is  discussed  in  the  AFWL  semiannual  report.^* 

A  different  and  very  simDle  technique  is  possible  using  the  triangular  interferometer,  and 
will  be  described  here.  The  relationship  between  the  differential  path  change  for  the  two  beams 
traversing  the  opposite  paths  in  fhe  triangular  interferometer  anii  ii>e  angle  between  ine  two 
beams  60  is  given  by1 


+  e2) 60 


(11-1) 


The  angle  60  is,  in  practice,  equal  to  the  pointing  angle.  (This  j  the  direction  the  incoming 
beam  to  the  interferometer  7.iakes  with  the  beam  that  would  traverse  the  interferometer  exactly 
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Fig.  II- 3.  Laboratory  boundary  layer  simulation:  (a)  ambient  air  temperature 
wind  (2  curves),  (b)  hot  air  temperature  wind  (6  curves). 
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symmetrically.)  In  operation,  the  shear  plate  is  continuously  rotated  at  a  fixed  speed,  and  the 
carrier  signal  frequency  produced  is  given  by  f  -  9A/9t  or,  from  the  above  expression, 


,  _  (H-  1 )  t  290 

Zp  9t 


Therefore,  keeping  00 /9t  constant,  depending  on  the  motor  speed,  the  carrier  signal  frequency 
is  propor  ional  to  the  pointing  angle  60.  If  the  interferometer  is  therefore  rigidly  fixed  in  a 
mount,  the  angle  of  arrival  of  a  wavefront  can  be  monitored  simply  by  measuring  the  frequency 
modulation  of  the  carrier  signal.  In  this  respect,  one  very  important  advantage  of  this  technique 
is  that  the  signal  amplitude  fluctuations  due  to  both  scintillation  and  to  phase  degradation  through 
the  atmospheric  path  do  not  influence  the  measurement  and  do  not  change  the  carrier  frequency. 

To  make  use  of  this  method  for  application  to  an  optical  tracker,  the  interferometer  pointing 
direction  can  be  controlled  by  a  servomechanism  o  reduce  the  carrier  signal  f  to  zero  and  keep 
it  there  so  that  the  pointing  angle  60  will  be  zero.  The  magnitude  of  the  minimum  angle  of  arri¬ 
val  (or  tracking  error)  that  can  be  anticipated  by  this  technique  can  be  estimated  by  using  Eq.  ( II- 1 ). 
The  n  inimum  angle-of-arrival  sensitivity  will  be  given  by 


where  A  can.  for  example,  be  taken  as  one-tenth  wave,  p  1.5,  t  =  2.5  cm,  and  0.  =0,  -  40”, 
min  12 

and  A  -  632.8  nm.  The  angular  sensitivity  is  thus 


2  sec  of  arc 


This  tolerance  is  calculated  for  the  case  where  the  interferometer  is  used  without  any  pre-optics 
(an  8-mm  input  beam).  Using  the  89-mm-diameter  Questar  telescope  produces  an  equivalent 
ten  times  magnification  in  the  angular  sensitivity,  whi'h  will  be  roughly  two-tenths  second  of  arc 
overall  sensitivity.  Longer  focal-length  optics  will  allow  correspondingly  larger  magnifications. 

It  is  clear  that  with  this  sensitivity  to  pointing,  rigid  mounting  of  this  equipment  is  mandatory. 

A  demonstration  of  this  technique  was  performed  in  the  laboratory  with  the  89-mm-diameter 
Questar  telescope  and  the  triangular  interferometer.  The  laser  source  (with  a  diverging  lens) 
was  located  about  14  meters  from  the  telescope.  A  small  0.95-cm-thick  optical  plate  was 
placed  close  to  and  in  front  of  the  source  on  a  rotatable  table.  One-degree  rotation  of  the  table 
produces  a  small  lateral  shift  of  the  source,  so  that  the  laser  beam  direction  as  seen  from  the 
telescope  changes  by  one  second  of  arc.  This  was  a  very  sensitive  way  of  changing  the  angular 
direction  of  the  laser  beam  into  the  Questar  telescope  and  interferometer  systems.  A  series  of 
interferometer  outputs  was  obtained,  shown  in  Fig.  11-4,  as  the  angular  direction  changed  through 
the  exact  "  on  axis"  position  and  beyond,  up  to  about  28  seconds  of  arc  off  the  axis.  The  carrier 
frequency  decreases  linearly  with  this  angle.  The  signal  frequency  modulation  is  zero  at  the 
on-axis  position.  The  width  of  each  picture  (each  double  curve)  in  Fig.  11-4  is  about  3  msec. 
Hence,  fast  response  in  an  FM  tracking  mode  is  quite  feasible.  There  is  some  dead  time  between 
curves  with  this  technique,  but  three  stacked  channels  in  one  interferometer  could  be  used  to 
overcome  this.  While  this  provides  tracking  along  one  axis,  a  similar  interferometer  stack  can 
be  used  to  provide  tracking  in  a  mutually  orthogonal  direction.  The  system  has  been  provisionally 
tested  in  the  laboratory  with  an  incoherent  light  source  (and  suitable  bandwidth  filters)  and  per¬ 
forms  satisfactorily  in  this  mode,  but  with  less  sensitivity.  At  a  later  date,  an  infrared  system 
will  be  tested. 
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Fig.  11-4.  Interferometer  output  for  different  angle-of-arrivals  (9'  (using  89-mm 
Questar  telescope). 
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